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Chapter 1 
Introduction 



The interplay between electron-electron interaction and electron-phonon coupling in 
strongly correlated fermionic systems has been intensively studied in the past and is 
still a major challenge for a realistic description and a better understanding of many 
materials. 

In the last decades, particular attention has been paid to one-dimensional fermionic 
systems, from both the theoretical and the experimental point of view. Starting from 
the features of one-dimensional electron systems, the initial goal was to understand 
correlation effects in higher dimensions, too. Furthermore, one-dimensional models are 
attractive from a theoretical point of view: there exists, in fact, a variety of methods 
that work exclusively in one dimension, sometimes even allowing for an exact solution. 
As a special feature of one dimension, the Fermi liquid theory breaks down and a new 
paradigm has to be introduced, namely the Luttinger liquid 13 • The first version 
of the generic model was proposed by Tomonaga, who showed that the excitations 
of a one-dimensional electron gas with linear dispersion are bosons, even though the 
elementary constituents are fermions. The excitations involve two particles, and the 
wave function of the two fermion states has bosonic properties. This model proved to 
be of fundamental importance for purely one- or quasi one-dimensional systems, e.g. 
quasi one-dimensional organic conductors [SI IH El El El or spin-Peierls compounds jSl E] 
where correlations are known to play an important role. In the last years, evidence 
for Luttinger liquid behavior has been found in semiconductor quantum wires fU] and 
carbon nanotubes [TT] . 

The concept of a Luttinger liquid is intimately connected with the bosonization 
technique whose origins date back to the seminal paper by Tomonaga £Q in 1950. 
During the following decades the method was worked out and successfully applied to 
one-dimensional electron and spin systems ^21 Ell El El- Despite its long history 
there are still some subtle points in the bosonization formalism which are not taken 
into consideration in the majority of the literature. One of these issues is the proper 
treatment of the so-called Klein factors which have to be introduced in order to pre- 
serve the anti-commutation property of the fermionic fields during the bosonization 
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procedure. Often, the existence of Klein factors has been ignored in the literature: 
this may be justified for infinite systems ^H], but in general they have to be treated 
carefully as pointed out, for example, in the context of impurity models and two-leg 
ladders [HJ EH 123 I2E] • The importance of Klein factors has also been emphasized 
by Schonhammer [2H1 I2ZI; who noted that the common practice of neglecting them 
may lead to erroneous results when nonlinear terms are considered. Such nonlinear 
terms arise in the presence of perturbations like impurity scattering or a modulation 
of the hopping. For the system with a perturbation an exact solution is known only 
in some special cases [HI]- In general one has to resort to approximative methods like 
renormalization group calculations [TZj. Another more intuitive method is the self- 
consistent harmonic approximation (SCHA) where the nonlinear terms are replaced by 
a harmonic potential with parameters to be determined self-consistently according to a 
variational principle for the energy or the free energy. The SCHA has been successfully 
applied to various nonlinear models [H d HH 121 121123 121 • 

In this work we handle the Klein factors in a systematic way, both in the thermody- 
namic limit and for finite systems. We develop an extension of the SCHA which treats 
the bosonic fields and the Klein factors on equal footing. As prototypical models we 
consider one-dimensional dimerized spinless fermions and Hubbard models with peri- 
odic modulation of the hopping (Peierls-Hubbard model) and of the chemical potential 
(ionic Hubbard model), respectively. Note that, throughout this work, we consider 
exclusively ground-state properties, i.e. the zero-temperature limit. 

We start with a short presentation of the bosonization method in chapter 121 where 
we focus in particular on the properties of the Klein factors. The chapter ends with a 
description of the Luttinger model, which, in certain cases, can be solved exactly. In 
the next two chapters we present a method for treating the Klein factors in bosonized 
Hamiltonians with nonlinear terms, and apply this method to spinless fermions in 
chapter El and to the Hubbard model in chapter 0] 

In chapter |5] we consider a model of one- dimensional interacting electrons coupled 
to a three-dimensional weakly correlated conduction band through a Coulomb inter- 
action. This model is motivated by materials like VO2, where at the Fermi energy 
one encounters both a band with one-dimensional dispersion and bands with three- 
dimensional character. We calculate the spin and the charge susceptibility in order to 
arrive at a description in terms of coupled one-dimensional chains with intra- and inter 
chain coupling in the spin and charge channel mediated through the 3d environment. 

In chapter El we summarize our results and give an outlook to some still open 
problems. Some technical details are described in six appendices. 



Chapter 2 
Bosonization 



2.1 Bosonization prerequisite 

Fermion systems in one dimension have features quite distinct from those in higher 
dimensions. In the early work of Mattis and Lieb j2H| and Bychkov et al. [SD] it has 
been shown that Landau type quasi-particles do not exist and the Fermi liquid theory 
breaks down in these systems. As a consequence a new concept has to be introduced 
to describe one-dimensional interacting electrons, known as the Tomonaga-Luttinger 
liquid. Corresponding models were first introduced by Tomonaga |]Q and Luttinger [2] 
and have been continuously improved in order to obtain an accurate characterization of 
real systems El CZl [T21 C31 H31 HSl The theory of the Luttinger liquid is closely 
related to the bosonization technique, which has been successfully applied to strongly 
correlated one-dimensional models. The basic idea behind bosonization is that particle- 
hole excitations have a bosonic character and are well defined quasi-particles at low 
energies. It consists essentially in a systematic mapping of a fermionic system (states, 
operators, Hamiltonians) into a bosonic one. It turns out that the bosonic language is 
often more suited for the understanding of the physics of the system, sometimes even 
allowing for an exact solution. 

A feature specific to one dimension and central to the whole development is the 
reduction of the Fermi sphere to two disconnected Fermi points ±k F . For energies 
close to the Fermi level, the dispersion relation can be linearized around the Fermi 
points. Following closely the notation of von Delft and Schoeller ^3], we consider a 
Hamiltonian with linear dispersion 

H o = ^ l e F-v F (k + k F )]c+ (7 c k(T + ^ {e F + v F (k-k F )]c+ a c kcT , (2.1) 

fc<0,cr fc>0,cr 

where e F is the Fermi energy and c\ a [c ku ) creates (annihilates) an electron with mo- 
mentum k and spin direction a =|, J.; we take h = 1 here and in the following. Since 
we are interested in low energy excitations, but not in the total ground state enery, we 
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vp{k + hp) + e F 




en(k) = v F {k - k F ) + e F 



Figure 2.1: Linearization of the spectrum around the Fermi points ±kp- The momen- 
tum domain of the left branch k G (—00, 0) and of the right branch k G (0, 00) was 
extended to k G (—00, 00). 

can extend the momentum domain of the first sum from k G (—00, 0) to k G (—00, 00), 
and the momentum domain of the second sum from k G (0, 00) to k G (—00, 00). In the 
following we will refer to these sums as left and right branches and we introduce the 
notion of left and right movers. As a result of the extension of the spectrum to minus 
infinity, the Hilbert space of the model is not the usual electron Hilbert space, but has 
been expanded to include a sea of unphysical states as well. This second, unphysical 
set of fermions requires high energies for their excitation and will not affect the low 
energy properties. In order to avoid the divergences associated with them, the normal 
ordering 



is introduced, where the reference state |0)o is the filled non-interacting Fermi sea. 
The physical fermion fields, defined as 



are separated in left and right movers with respect to the ±kp points as follows: 



:A: = A- o(0|A|0) 



(2.2) 




a k=— 



00 



(2.3) 



1 



00 



Ipphys 



VI 



a k>0 




(2.4) 
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For an interacting model, the separation into left and right moving fermions allows the 
sorting of various scattering processes according to their initial and final states [T7j 
("g-ology"). Using the notation tp kLa = c_ k _ kpa for the left movers and ^) kRa = c k+kpa 
for the right movers, the corresponding left and right fermion fields are defined by 

<Mx) = -= e ~ ikX ^L„ (2.5) 

k=—oo 



1 oo 

^Mx) = -Jf e< Nw, (2-6) 

k=—oo 

where again we extended the domain of momentum k e (—kp,oo) of each branch to 
k G (— oo, oo). These fermion fields satisfy the usual anti-commutation relations 

{^(x),^(x')} = 5(x-x')5 lY , (2.7) 

where 7 , 7 ' <= {R |, R j, L t, L j}. 

The density operator for particles with spin a is generally 

= J2 C t+ q a C kv (2-8) 

fe 

When the momentum q is small, it may be written as a sum of left and right movers 

Pa{q) ~ J] c fc+g- c ^ + Y C t+qa C ko 
k>0 k<0 

= Y ^k+qR^kRa + Y ^k-qLa^PkLa = pRa(q) + pLa(q)- (2.9) 



The left and right density operators PR a {q) and PLo(q) obey the following commutation 
relations: 



qL 

[pRa(q),PR.a'(-q')] = —^Saa'Sqq', (2.10) 



qL 

[pLa(q)ipLa>(-q')]= —^aa'^qq', (2.11) 



[PLa{q),PRa'{q')\ = 0. 



(2-12) 
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In evaluating the commutators (j2.10|) - (j2.12j) we use Eq. ()2.2|) for the density operators 
Pi(q) = '■ Pj(q) '■ ~~ o(0|p 7 (g)|0)o. The left and right charge and spin density operators 
are defined as 

pL/n{q) = pL/R]{q) + PL/Ri{q), (2.13) 

<TL/R(q) = PL/R]{q) - PL/Rl(q)- (2.14) 

In terms of left and right movers, the kinetic energy part of the Hamiltonian reads 

oo 

a k=— oo 

= V) / dx [■ ^RA x ) id xipRa(x) : + : ipt a (x)id x tp La (x) :] . (2.15) 

As a consequence of the linear dispersion, H and the density operators obey the 
following commutation relations: 

[H ,PRa(q)] = v F qp Ra (q), [H Q , p La (q)] = -v F qp La (q), (2.16) 

which allows expressing the kinetic energy part of the Hamiltonian as a quadratic form 
of charge and spin density operators. 

Using the density operators, one defines the following operators: 

b qLu = -^p La (q), b+ La = — PL{T (-q), (2.17) 



b qRa = -^ PRc7 (-q), b+ Ra = ^p Ra (q), (2.18) 

where n q = qL/2ir, and q > 0. These operators fulfill the bosonic commutation rules 

where 7 , 7 ' e {R T, R j, L T, L j}. 

Returning to real space, bosonic fields are defined as 

Vl«{x) = V -^e-^ a i/% La , (2.20) 
<PR*(x) = -J2 -^e iqx - aq/ %R«i (2-21) 

q>0 V U q 



2.1 Bosonization prerequisite 
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/ 












/ b qR(7 ~ 4>t-qRoi>kRa 



Figure 2.2: 6^ Lo . creates a particle-hole excitation on the left side, while b q R a annihilates 
a particle-hole excitation on the right side. 

where a — > is assumed. Formally, a~ l is a cutoff for the g-summation; physically, a 
corresponds to the lattice constant. It has been shown that the fermionic field operators 
are related to the bosonic fields via the following bosonization identity [T2*l IT^j : 

% La {x) = —F Lff e-*^ x) e- iVL ° ( - x) e-* riNL * x/L (2.22) 



i 



_FL a e~ ^" " e~ LaX/ , (2.23) 

^(z) = — ^F^e^^^+^^^^e 2 "^^. (2.24) 
V27ra 

The factor na in Eq. ()2.23|) comes from the commutator 

[<p+(x),(p y (x)]/2 = - In ^L/2na. 

The operators introduced in the bosonization identity count the number of 7- 
electrons relative to the reference state |0)o: 

00 00 

^7 = E : ^7^7 : = E Mr^*r - oW£lM>o] • (2-25) 



k=—oo k=— oo 



The _F+,F 7 operators, known as Klein factors or ladder operators, raise or lower the 
fermion numbers and assure that fermions of different species anti-commute. A detailed 
discussion of their properties is given in section 12.21 
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Later it will be convenient to introduce new charge and spin fields, <fi C:S and 6 C:S , 
which are linear combinations of left and right movers 

<Pc,s = TT/^^T ±( PH +<Pr\ ±<Pri + h.c), (2.26) 

9c,s = ^(mi^I-mTm + M, (2-27) 

where the index c stands for charge and the index s stands for spin. They satisfy the 
commutation relations 

M*)>My)] = Mz),9fy(y)] = 0, (2.28) 



Z7T ITT 

[4>/3(x), 9p(y)\ = -inQ(x - y)6pp> + —6 P p + — (x - y)8pp, a -»• 0, (2.29) 



%7T 

[(j)p{x), d y 0p,{y)\ = in5(x - y)5 pi3 > - — <W, a -> 0, (2.30) 

where 0,0' G {c, s}. In the thermodynamic limit, i.e. when the term \&pp> can be 
neglected, the fields 4>p{x) and d x 9p(x) are conjugate variables. 

We express now the electron density using the bosonization formalism. In real 
space, the density operator reads 

p„(x) = Y,J qx pM)~ z iqx P a(q)+Y, Yl e i{±2kF+q)x P*(±2k F + q), (2.31) 

q |g|<2fc F ± |<j|<2fc F 

which in terms of bosonic fields is 

Pa{x) « — — + — — + ^-d x [(f) La (x) + <j>Rj(x)\ 

L Z7T 



2na 



-2ik F x 2wi(N RlT +N L<T +l)x/L J^La+^Ra) „-ix(G-4k F ) 



_|_ ^ p+ p L ^ e -2ik F x e -2TTi{N RtT +N La +l)x/L Q -i{<j> La +<f> RtT ) _|_ ^ C M 32) 

Here, the term in the second line is important for commensurate band filling, when 
the reciprocal lattice vector is G ~ Akp. In the literature [20], ignoring both the Klein 
factors and the term in the second line, the density operator is simply given as 

P*( X ) = ^-d x [(j) LcT (x) + (j)Ra{x)\ + — COs[-2fc F X + 4>La( X ) + <f>Rcr(x)}- (2.33) 

Zix na 



2.2 Klein factors 
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F|0) = |-l)o F+\ - 2) = | - l)o F|3)o = |2) 

Figure 2.3: The action of F on the states |0)o and |3)o, and the action of F + on the 
state I — 2) . 



2.2 Klein factors 

The Klein operators F 1 and F^ lower or raise the total fermion number by one, which 
no combination of bosonic operators can ever do, and they assure also that fermions of 
different species anti-commute. In the thermodynamic limit and for gapless systems, 
for most practical purpose, e.g. the computation of correlation functions, there is no 
difference between states containing N 1 and N 1 ± 1 particles. However, when gaps 
open up, giving rise to finite correlation lengths, one must be careful when dealing 
with these operators [2U 123 HH] • 

For systems with only one type of fermions the Klein factors F 1 , simply modify 
the number of particles corresponding to the ground state on which the boson excita- 
tions are created, as depicted in Fig. 12.31 The Klein factors commute with the bosonic 
operators b gi , £>+ but not with the particle number operator N y : 

[F y , AT 7 ] = F y , [F+ N 7 ] = -F+ (2.34) 



If more than one species of fermions is considered, e.g. two in the case of spinless 
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fermions, the corresponding (Nl, iV^-particle ground states are tensor products of left 
and right movers: 

\N L ,N R ) =\N L ) ®\N R )o. (2.35) 

In order to preserve the fermionic character of the annihilation/creation operators, the 
Fr/Fr operators are defined such that they pick up a minus sign when they pass the 
left fermionic operators: 

F R \N L ,N R ) = (-1) Nl \N l ) ®\N r -1) 0j (2.36) 
F+\N L ,N R ) Q ee (-l)^\N L ) ®\N R + l) . (2.37) 

In general, for fermionic systems consisting of M species, one has to define a particular 
order for the species (Ni, Nm), such that 

17V,, N 2 , N M ) ee 1^)0 ® |iV 2 ) ® .... ® \N M ) (2.38) 

and the Klein factors pick up the total sign change from the preceding fermionic oper- 
ators: 

7-1 

F^\N X , N 2 , ...N M ) ee (-l)£i Na \N u N 2 , N y — 1, N M ) , (2.39) 

7-1 
V N <* 

F+\N h N 2 , ...N M ) ee (-l)-i 1^, iV 2 , iV 7 + 1, jV M ) . (2.40) 

For fermions with spin 1/2 (i.e. M = 4), the following order is defined: 

\N Lh N Lh N Rh N Rl } . (2.41) 

According to their definition, the Klein factors of different kind anti-commute. The 
commutation relations fulfilled by the Klein factors are 

fom F Y] = iKv F 7'] = [&<m *?] = *?] = °> ( 2 - 42 ) 

{F+ )J P y } = 2V, (2-43) 

{F 7 ,F y } = {F+,F+} = (if 7^7'), (2-44) 

[F 7 , 7V y ] = 5 7 yF 7 , (2.45) 



(2.46) 



2.3 Luttinger model 
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2.3 Luttinger model 

As a result of the approximations made, namely the spectrum linearization around 
the Fermi points and the introduction of a filled Fermi sea, some one-dimensional 
models become exactly solvable in certain cases, within the bosonization method. The 
Hamiltonian under consideration is 



H = H + H( 9l ) + H(g 2 ) + H(g 3 ) + H(g 4 ), (2.47) 



where H is the kinetic energy given in Eq. (j2.15|) and H(gi), H(g 4 ) are interaction 
terms, given by 



kpqcrcr' 

H(g 2 ) = jr ^ (92\\Saa> + g2±b-<Ta')^t Ra ll)p L(T 4p- q La>1pk- q R^ 
kpqaa' 

H(g 3 ) = ^7 (9n 5 ™' +93±$-aa>) 

kpqaa' 

X (^La^pLa'^q-pRa'^-k-q-(G-4k F )Ra 

+ ^ma^pRa'^-p-<l-(G-<lk F )La>4>-k+qLa), 

H{gi) = — Y (9i\\^aa' + ga.8-cr<T') 

kpqaa 1 

X (^kRa^pRa'^P+q^k-qRa + ^La^pLa'^p-oLa^k+qLa) ■ 



(2.48) 



For a graphical representation see Fig. 12.41 

The nomenclature in terms of gi is standard in the literature and is called "g-ology" 
|17j . The #2 and g 4 interactions are scattering processes with small momentum transfer, 
i.e. <C 2kp, while the g\ interaction is a scattering process with a large momentum 
transfer m 2kp. Umklapp scattering is relevant only for G ~ 4/cp, 6kp, where G is a 
reciprocal lattice vector, i.e. for commensurate filling. If the interactions are assumed 
to be instantaneous and local, the Pauli principle forbids the processes corresponding 
to g 3 |j and g 4 \\. In the absence of the g\±_ and g%±_ processes the bosonization of the 
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Figure 2.4: Graphic representation of the scattering processes. The solid and the 
dashed lines correspond to electrons belonging to the branches containing —hp and 
+kp, respectively. 



2.3 Luttinger model 
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terms ()2.15|) and (|2.48|) leads to the exactly solvable quadratic Hamiltonian 



H 



Luttinger 



E 

q>0 

E 

q>0 

E 

q>0 

E 

9>0 



( -rruF_ 9i\\ + g4± 
I L 2L 



-9i\\ + 92\\ + 92± 
2L 



( WVp_ g4|| ~ 94± 

1 L 2L 



(pL(q)pL.(-q) + pn(q)pR(-q)) 
(pL{q)pn{-q) + pR{q)ph{-q)) 

^L(q)(T L (-q) + <TR(q)(TR(-q)) 



-9i\\ + 92\\ ~ 92± 
2L 



(cr L (q)a R (-q) + a R (q)a L (-q)) 



+ "zero modes" (q — 0), 



(2.49) 



which in terms of dual charge and spin fields, Eqs. (j2.26)) - (J2.27)) . is known as Luttinger 
Hamiltonian 



H 



Luttinger 



E 

a=c,s 
7i 

AL 



E 



L dx 
2^ 

v, 



-(d x 4> a ) 2 + v a g a (d x 6 c 



N a + v a g a J a 



9a 



(2.50) 



The coefficients v a and g a are called Luttinger parameters, and they are related to the 
g's via [12] 



( 1 + J^L_^_) 2 _( 

\ 2-KVp 271V p J \ 



92\\ - g\\\ - g2± 

2ttv f 



1 , £4jj 5-4± 



2tivf 2-kvf 



ggjj ~ + ggl 
2^17^ 



27Ti; F + g4jj - gf 4 _L - g2|| + gi|| + g2± 
27TV F + g 4 || - g4± + #211 - glU - g2± ! 

1 2HV F + g4|| + g4± - g2|| + gl|| - g2L 
27TW F + # 4 || + g 4± + # 2 || - gl|| + g2_L ' 



(2.51) 
(2.52) 
(2.53) 
(2.54) 



The Luttinger Hamiltonian has the same type of low energy excitations as the harmonic 
chain, i.e. it corresponds to a free boson model (see appendix [BJ) . It contains two 
decoupled sectors corresponding to charge and spin excitations. The phenomenon of 
spin-charge separation is an important feature of fermionic systems in one dimension. 
For example, a "real" electron which is injected into an interacting system will decay 
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into its constituent elementary excitations which are charge and spin density modes 
that propagate at different velocities (v c and v s ). 

Besides the bosonic part, the Hamiltonian contains also a second part (given in the 
second line of (|2.50Jl ) which is a combination of N v operators, defined as follows: 



total charge: N c = N LA + N L>1 + N RA + N RA , (2.55) 

charge current: J c = N LA + N L>1 - N RA - N Rd , (2.56) 

total spin: N s = N LA - N L>1 + N RA - N R>1 , (2.57) 

spin current: J s = N LA - N LA - N RA + N Rtl ; (2.58) 

the fermionic ground state \N L ^, N L ^, N R ^, N R ^) can be redefined in terms of these 
new operators as 

\N Lh N Lb N Rh N Rl } = \N e , J c , N s: J s ) . (2.59) 



Chapter 3 
Spinless fermions 



In the Luttinger Hamiltonian, the Klein factors appear only in the combination F+F v = 
1 and hence they can be ignored. The Klein factors become important when nonlinear 
terms appear. In order to study their role in more detail, we consider in the following 
a system of interacting spinless fermions with nearest-neighbor interaction on a one- 
dimensional lattice with static dimerization. The static dimerization is introduced as a 
periodic lattice distortion, which leads to a modulation of the hopping amplitude. Both 
the nearest-neighbor interaction and the static dimerization introduce nonlinearities in 
the bosonized Hamiltonian, with non-trivial combinations of Klein factors F^Fr, which 
require a careful and rigorous treatment (321 • Using a variational ansatz in which the 
bosonic fields and the Klein factors are treated on equal footing, we calculate the 
energy gap and the Drude weight. To verify the accuracy of the method we compare 
our results with numerical and exact results which are available for certain values of 
the model parameters. 



3.1 Model and formalism 

The system of interacting spinless fermions with static dimerization is described by the 
following Hamiltonian: 

j 

+ V '52 n 3 n 3+l> (3.1) 

3 

where u is the dimerization parameter that leads to a periodic modulation of the 
hopping amplitude, and V is the strength of the nearest-neighbor interaction. In 
momentum space fermion operators are defined by 

cj= -= Ve ifcHj Cfc, with.Rj=ja, (3.2) 

V L / a k 
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Spinless fermions 



ID atomic chain 



Peierls distorsion 
* • # • 



Figure 3.1: Modulation of the hopping amplitude 



and correspondingly the Hamiltonian (J3.1|) reads 

H = e(k)clc k + 2itu ^ sin(/ca)c^c fc+7r / a 



k k 
q kp 



(3.3) 



where V(q) = Ve~ iga , e{k) = — 2tcos(ka), and a is the lattice spacing. 

In order to have a non-degenerate ground state we consider a half-filled system 
with an odd number of particles. Note that in this case, kp = ir/2a — n/L. Before 
bosonizing the Hamiltonian we have to choose a reference state with respect to which 
the boson excitations are created; we choose the Fermi sea of the half-filled system as 
this reference state, as illustrated in Fig. 13.21 In order to bosonize the Hamiltonian ()3.3)) 
we follow the notation introduced in chapter El The left and right spinless boson fields, 
and the bosonization identities are given by Eqs. ()2.20j) - (j2.23|) . where the spin index 
has to be dropped. Correspondingly, the charge and spin fields introduced as linear 
combinations of the left and right movers, given by Eqs. ()2.26)l - ()2.27j) . are replaced for 
spinless fermions by 



9{x) 



'■X<Pl-<Pr + <Pl-<Pr)- 



(3.4) 
(3.5) 



The bosonized Hamiltonian of this model reads 



H = H n + H ] 



Pcicrls 



Umklapp • 



(3.6) 
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Figure 3.2: Reference state for a system with 10 lattice sites, chosen as the filled Fermi 
sea for 5 particles. 
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+V / dx F+F+F L F L e~ Ai ^ x) e ~ ii7rxN/L + h.c. 
•/ o 



(3.7) 
(3.8) 
(3.9) 



where u = tu/(ira), V = V(2k F )/(2ir) 2 a and V(2k F ) = -Ve 2iwa ^ L « -1/. and are 
defined in Eqs. ()3.4|) - ()3.5j) and the cut-off parameter a is chosen equal to the lattice 
spacing. N and J are the charge and current operators defined as N = N F + Nr and 
J = Nl — Nr = 0, ±2, ±4, ... . The particle number N is a conserved quantity and is 
counted relative to the reference state, i.e. N = in the case of half filling. The g and 
v are the Luttinger liquid parameters, given in the weak coupling limit by 
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(3.10) 



where we have introduced the notation g 4 = V(0), gi = V(0) — V(2k F )- At half filling, 
for the model of spinless fermions with only nearest neighbor interaction, the Luttinger 
parameters can be exactly determined by comparison with the Bethe ansatz solution 



22 



Spinless fermions 



v/t 


-2 


-V2 


-1 





1 


i— 

V2 


2 


V 





tt/8 


7r/6 


tt/4 


tt/3 


3tt/8 


tt/2 


9 


oo 


2 


3/2 


1 


3/4 


2/3 


1/2 


v/ta 





0.94 


1.29 


2 


2.59 


2.82 


7T 



Table 3.1: Numerical examples for the relation between standard parameters of the 
spinless-fermion model at half filling, on the basis of the Bethe ansatz results in the 
regime —2t < V < 2t : V = — 2tcos(2r/), g = 7r/4?7, and v = irta sin(2r/) / (it — 2rj). 



|48j . The renormalized Fermi velocity v = irta sin(2r]) / (n — 2rj) and the Luttinger 
parameter g = 7r/4?7 are related with the interaction according to V = —2tcos(2rj). 

The Klein factors and the term proportional to the current operator irvgJ 2 / (2L) do 
not commute. We define a new operator A = F^Fl and its hermitian conjugate A + = 
F^Fr. The four-fermion terms arising from Umklapp scattering are F^F^FlFl = 
—A 2 and F^F^F^Fr = — (A + ) 2 . Since the Klein factors are unitary operators, it is 
easy to show that also the operator A is unitary and one concludes that its eigenvalues 
are pure phase factors, i.e. 

A\k) = e ik \k), A+\k) =e~ ik \k), (3.11) 

with < k < 2n. 

In the thermodynamic limit, L — ► oo, the term ~ J 2 can be neglected if we restrict 
the range of J values to a finite interval. We can choose a basis which diagonalizes the 
Klein factors and replace them by their eigenvalues. At the end one obtains 

tfpeieris = -uf dx itf k e" 2 ^ (x) Q~^ N I L + h.c. , (3.12) 
Jo 

iWiapp = +V [ dxe M e- 4! ^ e- 4l "^+h.c. (3.13) 
Jo 

As a result the Hamiltonian separates into different sectors of purely bosonic Hamilto- 
nians, labeled by k. 

As discussed e.g. by Schulz [T3], in the thermodynamic limit the Klein factors 
may be replaced by the Majorana ("real") fermion operators r\ w which assure the 
proper anti-commutation between the different species of fermions, but do not change 
the particle numbers of the ground state. The Majorana fermions satisfy the anti- 
commutation relation 



(3.14) 
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where a, (3 could be R, L and a for fermions with spin. A typical fermion interaction 
term, involving Klein factors, looks like 

ijj^ijjtipyips = FaFpFyFs x (boson operators). (3.15) 

Using Majorana operators a state with N particles is equivalent to a state with N ± 1 
particles and a fermion interaction term looks like 

^ipfip^ips = VaVpViVs x (boson operators), (3.16) 

with 

(VaVpVyVs) 2 = 1 VaVpV-fVs = ±1. (3.17) 

Using Eqs. ()3.14j) and (|3.17|) . one obtains only the eigenvalues r\ R r\ L = —r] L r] R = ±i, 
instead of e ±lk , i.e. the continuity is lost. Within the approximate approach using 
Majorana fermions the Hamiltonian (|3.fij) can be written as 

H = H + 2u dxcos20(x) -2V / dxcos40(x). (3.18) 
Jo Jo 



3.2 Phase diagram 

The zero-temperature phase diagram of spinless fermions with static dimerization has 
been investigated both numerically and analytically; numerically it is often more con- 
venient to study a dimerized spin 1/2 chain, which however can be mapped onto the 
model of interacting one- dimensional spinless fermions by the Jordan- Wigner transfor- 
mation [22] . In the V — u diagram one obtains different phases: a metallic one where 
the system behaves as a Luttinger liquid, and an insulating one where a charge gap 
exists. For u — > the critical point V/t = —a/2 that separates the two phases has 
been obtained in j^lj. For V/t £ (— v2, 2) the Umklapp scattering process is irrelevant 
and it is renormalized to zero, while for V/t > 2 this process becomes relevant and 
results in the opening of a correlation gap. For u^O, according to perturbative renor- 
malization group calculations (22] and confirmed numerically using the density matrix 
renormalization group (DMRG) [2H]; the phase boundary in the V — u plane shifts to 
more negative values of V with increasing u, as can be seen in Fig. 13.31 



3.3 Self-consistent harmonic approximation 

In the following we consider the bosonized form of the model given in Eq. (|3.6J) . To 
discuss the gap formation, we use the self-consistent harmonic approximation (SCHA) 
which has originally been introduced for the sine-Gordon model jTHj. The idea is to 
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Figure 3.3: Phase diagram for dimerized spinless fermions as obtained using the DMRG 
|23j . The metallic phase of the system is denoted by o. As demonstrated in [23] the 
metallic region increases when including not only a periodic modulation of the hopping 
parameter, but also in the interaction (+)• The region with a gap in the excitation 
spectrum is marked by □. 

construct a trial Hamiltonian H tr where the nonlinear terms ~ e ±2tr t> [ n Eq. ()3.8|) are 
replaced by a quadratic form with a parameter A to be determined self-consistently 
according to a variational principle for the energy. 

As we discussed in the preceding section, in the thermodynamic limit we can replace 
the Klein factors by their eigenvalues Eq. (|3.11|) . where we choose the phase factor 
k = 7r/2. The resulting Hamiltonian is equivalent to the Hamiltonian obtained within 
the Majorana fermion approach and is given in Eq. (|3.18|) . 

We introduce the trial Hamiltonian 

K = f L ^{ V -{d x d ) ) 2 + vg{d x ef + —^{x)\. (3.19) 
Jo ^ 19 vg J 

Since we are only interested in properties at zero temperature we use the normalized 
ground state \i[}q) of H tT as a variational wave function for the Hamiltonian ()3.6|) . This 
provides us with an upper bound E for the ground state energy E of ()3.6|) due to the 
inequality 

E < E = {MH\^)- (3.20) 
The variational energy is given by 

E = (H )ti + (#Peierls)tr, (3-21) 

E EL A 2 . „. 
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where E tT is the ground state energy of H£. Minimizing E with respect to the varia- 
tional parameter A yields the gap equation 

A 2 = Snvgu e -2 ^* 1 . (3.23) 

Since H£ is bilinear in the field operators it is straightforward to calculate the equal 
coordinate correlation function of the phase field 4>(x) entering Eq. (|3.22|) 

where the k- values in the sum are multiples of 2ir/L. Replacing the sum in Eq. ()3.24|) 
by an integral we obtain (for A <C Ao) 

(0 2 )t r = fln^, (3.25) 

where Ao = Iv^T 1 ja and 7 = 0.5772 is Euler's constant. Inserting Eq. ()3.25|) into Eq. 
flH2BD yields 

A ( u\ 1,{2 ~ 9) 

#"=(-) > (3-26) 

A \u J 

with Uq = e~ 7 A /4gl For g > 2 the right hand side of Eq. (|3.26|) diverges as u — > 0, 
and the equation has only the trivial solution A = 0, i.e. the line g = 2 marks the 
transition from a gapless to a gapped phase in the g — u plane. In the spinless fermion 
model this corresponds to the line V/t = — v^2 in the V — u phase diagram. The 
value g c = 2 obtained within the variational approach is in accordance with the exact 
result |3*lj . However, according to renormalization group (SHj and DMRG calculations 
|23j . the phase boundary in the g — u plane shifts to larger values of g with increasing 
u, while the SCHA gives a vertical line. On the other hand, the exponent 1/(2 — g) 
characterizing the opening of the gap is exact [T§] . 

For non-interacting fermions V — 0, the Luttinger parameter is g — 1 and corre- 
spondingly A a it, in agreement with the exact energy gap. The prefactor however 
depends on the cut-off procedure used. From Eq. (|3.26J) we find A = 4tMe 7 , where the 
factor e 7 has its origins in the soft cut-off used in Eq. (|3.24j) . If we replace in Eq. ()3.24|) 
the soft cut-off by a hard cut-off the result is A = Atu/n. 

At V/t = 2, g=l/2, where Umklapp scattering becomes relevant and leads to the 
opening of a correlation gap, Eq. (J3.26|) yields A oc u 2 ^ 3 which agrees up to a logarithmic 
correction with the exact result A ex oc w 2 / 3 / yj\ h\u\ 
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3.4 Finite systems 

In the following we consider a finite system of length L. In this case it is not possible to 
replace the Klein factors by their eigenvalues, since the term proportional to the current 
operator J 2 in the Luttinger Hamiltonian (|3.7|) does not commute with the F's. The 
trial Hamiltonian is now chosen as a sum of two commuting parts H tT = + H^, 
with given in Eq. ()3.19|) and the "Klein Hamiltonian" 



H* = -iB{FiF R - F+F L ) + ^ J 2 , (3.27) 

where A and B are variational parameters to be determined self-consistently. is 
of the form of a tight-binding Hamiltonian for a particle moving on a one-dimensional 
lattice in a harmonic potential. A class of similar Hamiltonians has been studied in [2*7j . 
In this case we decouple the Klein factors from the bosonic fields and the minimum 
condition for the variational energy is equivalent to the following substitution in the 
Peierls dimerization term 

F+F R e 2l ^ x) -> (F+F R )e 2l ^ x) + F+F R (e 2 ^ (x) >. (3.28) 

As a result, instead of replacing the products of Klein factors by their eigenvalues as 
in the thermodynamic limit we now have to replace them by their expectation values 
with respect to the ground state of H^. According to the Feynman-Hellmann theorem 
(HZ! EE], the expectation value (F£F R ) can be expressed in terms of the derivative of 
the ground state energy E of the Klein Hamiltonian 

(F£F R ) = z||. (3.29) 

The variational energy obtained from H tr has to be minimized with respect to the 
parameters A and B and we obtain the gap equations 

A 2 = -ATxvguE'(B)e- 2 ^\ 

(3.30) 

B = uLe~ 2 ^, 

where E (B) is the ground state energy of and E' (B) = dE (B)/dB. In order to 
calculate the ground state energy of it is convenient to switch to the momentum 
representation, J — > — id/dp j, where the Hamiltonian reads 

H * = ^-aW J +Wc0s{2pj) > (3 ' 31) 

with periodic boundary conditions for the wave function, ^>(pj + vr) = ty(pj). At this 
point, we note a close analogy with the quantum theory of Josephson junctions, see 
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Figure 3.4: Energy gap A(L) (in units of t) as function of the dimerization parameter u 
for L = 100, 200, 500, 1000, 2000 (from top to bottom). The dashed line is the analytic 
result of Eq. (|3.26|) valid for L —> oo and A -C A . 

js23 and appendix O The corresponding Schrodinger equation is the Mathieu equation. 
Accordingly, the ground state energy of Hf r has the following asymptotic behavior |4Ti] : 



E (B) « <^ WVg /^— (3-32) 



Combining the two equations (|3.30|) one obtains 

A 2 L 2 = —4irvgE' (B)BL 1 (3.33) 

which can be used to relate the two limiting cases of Eq. ()3.32|) with two different 
physical situations: The condition of large system size L ^> v/B is equivalent with 
A ^> v/L and therefore the results for (0 2 )t r and A (it) are the same as for the infinite 
system. On the other hand, the condition of small systems, L v/B implies that 
A v/L. In this case the size dependent energy gap A(L) is purely due to the finite 
system size, i.e. it is of order v/L. Thus the crossover from a finite-size gap to a true 
dimerization gap coincides with the crossover between the regions where Klein factors 
are relevant or can be ignored. 

In general, the set of equations (|3.30J) can only be solved numerically. Fig. 13 .41 shows 
the size-dependent energy gap A(L) = a/ (2ttv/L) 2 + A 2 as function of u for several 
values of the system size L; notice that periodic boundary conditions are imposed, 
where the smallest wave number is Itx/L. The Luttinger parameter is (7 = 3/4 which 
corresponds to V/t = 1. 
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3.5 Drude weight 

The Drude weight or charge stiffness is a physical quantity used to characterize charge 
transport at zero frequency. At T = the real part of the electrical conductivity at 
frequency u is of the form ct{uj) = 2ttD5{uo) + <j reg (uu) with a reg (uj) — > for uo — > in 
a system without impurities. Therefore D = characterizes an insulator while D > 
describes an (ideal) conductor j^TJ. The Drude weight is not a true order parameter 
related to a symmetry breaking, however, it can be used to distinguish the insulating 
from the metallic phase. The simplest way to calculate D is due to Kohn (3^ who has 
shown that the Drude weight can be expressed as 



D _L d 2 E(tp) 



2 dtp 2 



v=0 2 V 27T J d& 



(3.34) 



$=0 



where E(tp) is the ground state energy of a ring of circumference L which is threaded 
by the magnetic flux $, where tp = 27r$/$o and $o = h/e. In the fermionic model 
(|3.1j) the hopping parameter t picks up a phase factor exp(±i(p/L) when a flux $ is 
applied. Alternatively, within a proper gauge transformation, the parameter <p can also 
be associated with a change of boundary conditions, i.e. tp — corresponds to periodic 
and (/? = 7T to anti-periodic boundary conditions. The Drude weight is directly related 
to the persistent current /($) which flows in a ring that is penetrated by a magnetic 

flux muni 



3.5.1 Drude weight within bosonization 

We now turn to the calculation of the Drude weight D within the bosonization for- 
malism. In the presence of the magnetic flux the only modification in the bosonized 
Hamiltonian of the Luttinger model appears in the current operator J which has to be 
replaced by J + - jH] and the Drude weight is given by the product of Luttinger pa- 
rameters Do = vg/2n. Correspondingly, for the system with dimerization, we modify 
the -B-dependent part of the the trial Hamiltonian and write 

= -iB{FiF R - F+F L ) + ^(J + ^) 2 - (3-35) 

Applying the same procedure as in the case tp — yields a variational estimate E for 
the ground state energy which now depends on the flux tp, i.e. E = E(A, B, tp) where 
A and B, obtained from the gap equations ()3.30|) . are also functions of <p. On first 
sight it might seem impossible to calculate the Drude weight using Eq. ()3.34|) since 
there is no analytical solution of the gap equations ()3.30|) even for tp = 0. However, a 
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closer look reveals a great deal of simplification. From the second derivative 



d 2 E 

dtp 2 



d 2 E /dA 
dE 8 2 B 



d 2 E /dB 



dE d 2 A 



dB 2 \dipj dA dip 2 



d 2 E 



one retains only 



d 2 E 



dp 2 



dB dip 2 dip 2 



d 2 E(A,B,<p) 



+ mixed terms 



(3.36) 



ip=0 



dip 2 



(3.37) 



ip=0 



dE 
dB 



due to the minimum condition of the energy and ii) c 



dip 

Since the 



since i) ^ 

^ = at ip = due to symmetry (A and B are even functions of cp) 
trial Hamiltonian H tT = H£ + (<p) consists of two commuting parts we may write 
E(A,B,ip) = E(A) + E (B,ip) where E(A) depends only on A and E (B,ip) is the 
ground state energy of H^(p>), i.e. we obtain the simple result 



D 



L d 2 E (B,ip) 
2 



dip 2 



(3.38) 



<p=0 



where B is given by Eq. ()3.30|) . 

To proceed we again represent H^(ip) in terms of the Mathieu equation (|3.31j) where 
now the boundary conditions are ^(pj + it) = e lip '^{pj). It is now straightforward to 
calculate the Drude weight for a finite system of size L in the gapped phase g < 2. In 
the "finite size gap" region (LA t>), the confining potential of the trial Hamiltonian 
(|3.35|) is the dominant term and the hopping can be regarded as a small perturbation. 
In second order perturbation theory we obtain 



D = D (l 



(3.39) 



where q = 2LB/nvg oc uL 2 ~ 9 and Dq = vg/2ir is the Drude weight of an unperturbed 
Luttinger liquid. In the limit q — > which is equivalent to uL 2 ~ 9 — > the Drude weight 
is equal to the value D . 

In the opposite limit (LA ^> v ) which corresponds to a Mathieu equation with a large 
cosine potential the variation of the ground state energy with change of boundary 
conditions is exponentially small (as expected from the WKB approximation). A more 
careful treatment (see appendix |Dj) yields the exponential dependence 



D pa D n — 



4 / LA\ 



7T 



3/2 



vg 



J 



-2LA/-rrvg 



(3.40) 
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where A and u are related via Eq. ([3.26)1 . 

For g > 2 the system remains in the Luttinger liquid phase for small dimerization 
u, however with renormalized Luttinger parameters. These renormalizations are not 
obtained within the SCHA and consequently the Drude weight is not reduced in this 
region. 

3.5.2 Drude weight for free spinless fermions 

As a test case that can be solved exactly, we consider the Hamiltonian (|3.1)1 in the 
limit V = 0. In this case we can calculate the ground state energy analytically. Our 
aim is to compare the exact result for the Drude weight with the result obtain within 
the SCHA in order to assess the accuracy of the method. 
When a flux $ is applied, the system is described by the Hamiltonian 

H = -t^(l + n(-l) J ')(c+c j+1 e^ /L + h.c), (3.41) 

3 

which in momentum representation reads 

H = -2t s T'(c + c+ )( C ° sk ~ iusink \( °k ) (3 42 ) 

Z_s v k ' k+nJ y iusink —cosk J \ Cfe +7r J v 1 

The prime on the sum indicates the reduction of the Brillouin zone from [—-,-] to 
[— The corresponding dispersion ei t 2(k) = ±2tV cos 2 k + u 2 sin 2 k is shown in 

Fig. 13.51 The ground state energy is 

E(<p) = -2t^2' Vcos 2 k + u 2 sin 2 k, (3.43) 

k 

where for a system of length L, k = (2irn + tp)/L, and the lattice constant is set to 
one. In the following we assume a dependence of the ground state energy on <p of the 
form 

E(<p) = E + E 1 coa(<p), (3.44) 

The second term, in general, is an infinite sum of harmonic terms, but for a ring with a 
Peierls distorsion, threated by a magnetix flux, at T = the main contribution is given 
by the first harmonic term [33] . Under this assumption the Drude weight is given by 

D = = j[E(n) - E(0)}. (3.45) 

We calculate now the difference AE = E(tv) — E(0) which is needed in the Drude 
weight evaluation (j3.45[) . We introduce z m = e lkm and the function f(z) = l/(z L — e lip ) 
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Figure 3.5: Dispersion of free fermions with dimerization u. 



which has single poles at z m with residua Kesf(z) 



2 e l(p /L. Expressing Eq. 



()3.43|) in terms of z m we may replace the sum by a contour integral and obtain 
tLe^Vl - u 2 



E{<p) 



dz \/ z 2 + z~ 2 
2ni z(z L - e-w) 



27 



(3.46) 



where 7 



u 2 )/{l 



u 



and C is a contour that encloses the singularities of f(z) 



We choose C to be composed of two circles with radii slightly larger or smaller than one, 
respectively. Substituting z — > l/z the integral along the inner circle can be mapped 
onto the integral along the outer circle. The square root in the numerator has branch 
cuts along the imaginary axis in the interval [—2/1,2/1] and for \y\ > yo where ±ij/o 
and ±m/i are the zeroes of the function under the square root. Now we deform the 
integration contour along the branch cut from ±ij/o to ±zoo, as can be seen in Fig. 13.61 
and obtain for the energy difference AE = E(n) — E(0): 



AE 



u 



2 r°° 



7T 



dy y / y 2 + y~ 2 - 27 



y y L -y L 



(3.47) 



Substituting y = e x yields 
AE 



AtL^T- 



7T 



, J 2 cosh(2a;) - 2 7 

dx ; -. , 



(3.48) 



with 2xq = Arcosh7 = ln((l + u)/(l — u)). For Lxq 3> 1 the integral is rapidly cut 
off by the exponential and we may expand the square root around x = xq. The Drude 
weight is then 



D 



LAE 



2uL 



7T 



2 1 — u 



1 



3(1 



u 



8Lu 



(3.49) 
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Figure 3.6: The contour C is deformed into the contour C along the branch cuts between 
±iyo and ±ioo. 

Expanding the logarithm for ti<l, and inserting A cx = 4tu and vp = 2t yields 



In the special case of free fermions, the Luttinger parameters are g = 1 and v = vp = 2t. 
If we insert these values in Eq. (|3.4()jl we can compare the exact result, Eq. (|3.5()j) . with 
the Drude weight obtained within the SCHA, Eq. ([3.40)1 . We see that the exponential 
behavior D ~ exp(— const • LA/v) is reproduced, however, with a different constant 
2/ 7i instead of 1/2. This slight difference is not surprising since the energy scale A , 
see Eq. (I3.26|) . is known to depend on the cutoff procedure within SCHA. 




(3.50) 



Chapter 4 
Hubbard model 



In order to analyze further the importance of Klein factors in bosonized Hamiltoni- 
ans with several nonlinear perturbations, we consider in the following a system of 
interacting spin 1/2 fermions on a one-dimensional lattice. As in the case of spinless 
fermions, we demonstrate how to handle the Klein factors in a systematic way, both in 
the thermodynamic limit and for finite systems. Within the self-consistent harmonic 
approximation the decoupling of the Klein factors from the bosonic fields results in 
a two-dimensional tight binding Hamiltonian with a confining potential. While for 
spinless fermions, the Klein Hamiltonian can be mapped onto a Mathieu equation 
which can be solved analytically in certain limits, for fermions with spin the prob- 
lem becomes more complicated. As an application we consider two different models, 
the Peierls-Hubbard model and the ionic Hubbard model. The first model contains a 
static alternating Peierls distorsion u of the lattice, i.e. a modulation of the hopping 
amplitudes in the kinetic term. In the second model, an alternating on-site energy 
modulation is introduced, i.e. a staggered potential. 

4.1 Hubbard model 

The Hubbard model is described by the Hamiltonian 



where cf a creates an electron with spin direction a =f, | and t is the hopping para- 
meter. U is the interaction energy of two electrons at the same site, (U > 0). 
In momentum space the Hamiltonian reads 






(4.2) 



34 



Hubbard model 




Figure 4.1: Schematic representation of backward scattering. Solid (dashed) lines refer 
to electrons close to the left (right) Fermi points. gi± couples electrons with opposite 
spin. 



where €k = —It cos(ka) and N is the number of lattice sites. Applying the bosonization 
procedure, described in chapter |21 one obtains 

H = H + H 1 + H 2 , (4.3) 

where H is the Luttinger Hamiltonian, 

-L 



Jo 2vr t g a J 4L ^ \ g a 



(4.4) 



and Hi and H 2 are the backscattering and Umklapp contribution, respectively. Back- 
ward and Umklapp scattering, which represent scatterings with large momentum trans- 
fer, introduce nonlinear terms in the spin and the charge sector, respectively. As a 
consequence the spin and the charge currents are no longer conserved. 

The Hamiltonian resulting from the backward scattering corresponds to an interac- 
tion between left and right Fermi points with a momentum transfer q ~ 2kp, see Fig. 
14.11 In bosonized form the backward scattering term reads 

Hi = U ( L dx iF+F+FR^e-^^e-^* 
Jo L 

+ F+F+ i F Li F xr e i! * N '*e 1 * /s +'], (4.5) 

with U = UL/[(27ra) 2 N]. 

The Umklapp term H 2 is only important for commensurate band fillings. This 
process comes about because the state at Uf is equivalent to the state kp + nG, where 
n is an integer and G is a reciprocal lattice vector, G = 27r/a, with a being the 
lattice spacing. For the half-filled band case, kp = vr/2a, and therefore the scatterings 
kF — > &f + Ikp and kp — > —kp are equivalent. The Umklapp process is characterized 
by scattering of two particles in the same direction across the Fermi surface with a 
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Figure 4.2: Schematic representation of Umklapp scattering. Solid (dashed) lines refer 
to electrons close to left (right) Fermi points. g%± couples electrons with opposite spin. 



momentum transfer q ~ 4k p, see Fig. 14.21 
The bosonized form of the Umklapp term is 



-L 

Ho = U I dx 



[ dx \F^F^F Ll F L ^ N ^-^ 
Jo 1 



+ F^F Rl F m e^ N ^ 2 ^\ . (4.6) 

The Hubbard model is a particular case of the "g-ology" model given by Eq. ()2.4|) . for 
which the coupling parameters are gi± = U and gq = 0. From Eqs. (j2.f)lj) - (J2.54j) the 
corresponding Luttinger parameters are obtained as 



U U 

V F \1 , V C = V F \1 



7TVp V KVp 



-— \A + — 

nvp y nvp 

where vf = 2tasin(/c^a) is the Fermi velocity. The ground state phase diagram of 
the "g-ology" model has been investigated in detail using renormalization group tech- 
niques [TJj. Schulz ^2] noticed that in the small U perturbative regime, for all filling 
factors the backward scattering renormalizes to the weak-coupling fixed point gl ± = 
and the g s parameter scales to g* — 1. In lowest order of the renormalization group 
equations the charge parameter g c remains unrenormalized and its perturbative value 
is g c = 1 — U /2ttvf 1 . Schulz also pointed out that there is no phase transition in the 
Hubbard model between weak and strong coupling, and hence the small and large U 
limits belong to the same phase. As a consequence, the equality g s — 1 remains valid 
for all repulsive interactions. 



1 There is a misprint in where the perturbative value of the charge parameter is given by 
g c = 1 - U/nv F . 
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In the strong coupling limit (U ^> t), there are no low-energy charge excitations involv- 
ing doubly occupied sites. In a restricted Hilbert space containing only singly occupied 
sites, in second order mt/U the Hubbard Hamiltonian can be mapped onto the t — J 
model, where the spins are coupled through an effective antiferromagnetic exchange 
integral ,7 = At 2 IU P] 

H t -j = -*y~] [(1 - n j _ .)c+.c j+ i CT (l + h.c.]+J ^ 

j i 

At half filling the t — J model reduces to the antiferromagnetic Heisenberg model 

H = J^S r S j+1 , (4.8) 

i 

which can be transformed to the spinless-fermion Hamiltonian with nearest neighbor 
interaction using the Jordan- Wigner transformation discussed in the previous 
chapter, for the special case V — 2t, g — 1/2. 

The exact solution of the one-dimensional Hubbard model via Bethe ansatz has been 
found by Lieb and Wu |lHj. However, the Bethe ansatz wave function is in many cases 
too complicated for the evaluation of expectation values or correlation functions. Still, 
the Bethe ansatz is a convenient starting point to determine the Luttinger parameters 
for arbitrary band filling and interaction strength |13j . 



1 

T 



(4.7) 



4.2 Peierls-Hubbard model 

The coupling of a one-dimensional metal to an elastic lattice results in an instabil- 
ity towards a lattice distortion known as Peierls instability [201 • A lattice distorsion 
decreases the electronic energy which overweights an increase in lattice energy. Due 
to an opening of a band gap at the Fermi energy the system becomes an insulator. 
This phenomenon has been observed in various quasi one- dimensional materials, e.g. 
conjugated polymers [7j like polyacetylene, charge-transfer salts [SJ and more recently 
in the one- dimensional spin systems CuGeOs [H] and NaX^Os 0. 

From a theoretical point of view, a first step towards a quantitative description of 
the Peierls transition has been made by Su, Schrieffer and Heeger j^2]- However, in 
their model no electron-electron interactions are taken into account. 

The effects of electronic correlations on the Peierls transition have been studied 
using a variety of methods including variational wave functions Hartree Fock 

quantum Monte Carlo [SH], numerical diagonalization of small systems |5T| . 
bosonization [Hill EI], and incremental expansion 155]. 
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In the following we neglect the lattice dynamics and consider a model with static 
lattice distortion given by the Hamiltonian 

H = -t£)(l + (-l)%)(c+ c j+la + c+ +la c ja ) + Uj2n^n n , (4.9) 

J> 3 

which is the usual Hubbard Hamiltonian with an additional periodic modulation of the 
hopping described by the dimerization parameter u. For the case of a half-filled band, 
which we consider here, the lattice modulation is alternating. In momentum space, the 
Peierls contribution reads 

tfpeicrls = tU ]T(e^C+ C fe+mT + h.C.) (4.10) 
k,(T 

and in bosonized form 



f dxF+F Ll e- i ^ N °- N ^ x e- i ^-^ + h.c.. (4.11) 
na J 

At half filling, N c = N s = 0, and using the notation u = tu/ira we obtain 

^Peicris = u dx {iF^F L ^~ l ^ + ^ + iF+ F^e"*^"^ + h.c.}. (4.12) 

./O 

In the following we focus on the role of the Klein factors. Due to the conservation of 
charge and spin all combinations of Klein factors appearing in the Hamiltonian (|4.12jl 
can be expressed in terms of the operators Af = F^Fl] and A± = F^F^ plus their 
hermitian conjugates. In particular, the four-fermion terms arising from Umklapp and 
backscattering read 

F ri F ri F liF L1 = F+F LT F+F Li = A^A b (4.13) 
F+F+F Rl Fv = F+F^F^ = A T A+. (4.14) 

Since the Klein factors are unitary, F+F a = F a F+ = 1, it is easy to show that 

[A h A+) = [A l ,A+] = 0, (4.15) 

and we may choose a basis where A a and A+ are both diagonal. From Aj!~ A|- = Al = 
1 one concludes that the eigenvalues of A a are pure phase factors, i.e. 

A^k x ) = e*^,*;), A+l^ki) = e-^l^,^), (4.16) 
A^k,) = e*^,*;), A+lfet.fc;) = e-^|fc T ,^), (4.17) 
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with < k a < 2ix. Although the terms ~ J% a appearing in ()4.4|) do not commute 
with the Klein factors it seems reasonable to neglect them in the thermodynamic 
limit L — > oo. We may thus replace the Klein factors in H\,H2 and -£/p e i er is by their 
eigenvalues, and obtain 



Hi = U dx { e i(h-h) e -xV2<f> s + h c ^ 
Jo 



(4.18) 



r L 

H 2 = U dx {e^t+^e-^v^c + h c ^ ( 419 ) 
Jo 



^Peierls = U \ dx {ie^e'^^ + feik^-iyftfc-M + h c J ( 4 _2 ) 

Jo 

As a result the Hamiltonian of the dimerized Hubbard model separates into different 
sectors of purely bosonic Hamiltonians which are labeled by and k^. In an approach 
where the Klein factors are replaced by Majorana fermions JH] one obtains only the 
eigenvalues ±z for the two-fermion terms, and ±1 for the four-fermion terms, i.e. the 
continuous symmetry is lost. 

4.2.1 Self-consistent harmonic approximation 

In order to study the opening of charge and spin gaps in the dimerized Hubbard model 
we use the self-consistent harmonic approximation (SCHA) in which the exponentials 
of field operators appearing in ()4.19j) - (j4.2()j) are replaced by quadratic forms. In the 
thermodynamic limit we replace the Klein factors by their eigenvalues and introduce 
the trial Hamiltonian 

t ^{-{dA a f + v a g a {d x e a ) 2 + ^<A, (4.21) 
^ s Jo 2tt [g a v a g a J 

which provides us with a variational estimate for the ground state energy 

E = (^o)tr + (#l)tr + (H 2 )ti + (#Peierls)tr, (4.22) 

j = + (4-23) 

a=c,s J 

where E tr is the ground state energy of H tI , and 

e(k h fcj.) = 2B c cos(A; t + fy) + 2B s cos(k ] - k±) - 25 cs (sinA; T + sinfy), (4.24) 
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range 


i 


; 

h 


( TD TD TD \ 


< B cs < 2B S 


arcsin jff- 


7r — arcsin jff- 


-2B C -2B s -§f 


2B S < B cs 


7T 

2 


7T 

2 


— 4-Bcs — 2£> c + 2£> s 



Table 4.1: Minimum eo of e(fe-j-, /cj) = 2B C cos(/c| + fcj) + 2B S cos(fc| — fcj) — 2i? cs (sin + 
sinfcj) which is used in the gap equations (J4.28)) and (|4.29j) . The first line is only 
relevant for u = (Hubbard model without dimerization) whereas the second line is 
relevant for u > 0. 



with 



B c = Ue~^^, (4.25) 
B s = f/ e - 4( ^ >tr , (4.26) 
B cs = ue'^e'^. (4.27) 



-AB C — §- - B cs —^, (4.28) 



Minimizing the variational ground state energy with respect to A c and A s yields the 
gap equations 

2nv c g c c dB c cs dB c 

= -4B.P- - B^pL, (4.29) 
27rv s g s dB s dB cs 

where Cq is the minimum of e(k^, k^) with respect to k^ and k^ (see Table l4~Tf . 

The first minimum is only relevant for u = where B cs = 0, i.e. the first line of 
Table l4~T1 with deo/dB c = —2, deo/dB s = —2 corresponds to the Hubbard model 
without dimerization. In this case the equations for the charge and spin parameters 
are 

A2 

8B C , (4.30) 
8B S . (4.31) 



2irv c g c 
2ixv s g s 



The numerical solution of these equations is displayed in Fig. 14. 3[ where here and in 
all the following figures both A c and A s are given in units of t. 

Since H tT is quadratic in the bosonic fields it is straightforward to calculate (4>l)tr 
and (</>s)tr analytically. In the thermodynamic limit, where the sums are replaced by 
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integrals one obtains 

(<f>X = fln^, (4.32) 
= fln^, (4.33) 

where A 0c ,os are cutoff-dependent energy scales of the order of the bandwidth, and 
A CiS < A OCj os is assumed. The solution of the gap equations ()4.3()j) and (|4.31j) is 

/ x V(2-2g c ) 

At = Ac0 , (4.34) 

A s = 0. (4.35) 

Using the perturbative result g c = 1 — U/2ttvf and vp = 2t, we obtain 

Acoce- 27 ^*/^ (4.36) 

in the small U limit. This should be compared with the exact solution from Bethe 
ansatz where for small U the charge gap is given by 

A B.a. = 8U e -2nt/U^ (4 37) 

Tit 

The SCHA result deviates from the exact solution by the factor \n(t/U) in the exponent. 
We show the charge gap obtained within the SCHA in comparison with the exact 
solution from the Bethe ansatz in Fig. 14.41 

For non-zero dimerization u > a solution of the gap equations exists only for 
B cs > 2B S , i.e. the second line of Table I4~T1 has to be used, with de /dB c = —2, 
de /dB s = 2 and de /dB cs = —4. The corresponding gap equations are 

A 2 

c — = AB CS + 8B C = Au e-^-W* + 8U e " 4< ^ )tr , (4.38) 

2nv c g c 



'— = AB CS - 8B S = Au e-^-W* - 8U e - 4( ^K (4.39) 

2nv s g s 

The numerical solution of these equations is represented in Fig. 14.51 for U/t = 1. For 
very small dimerization u the charge gap A c approaches a constant while the spin gap 
goes to zero, as expected for the Hubbard model without dimerization. On the other 
hand, for strong dimerization A c and A s approach each other as expected for a band 
insulator without correlations. The transition between the two regimes is smooth and 
continuous. 
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Figure 4.3: The charge gap A c (in units of t), obtained within SCHA as function of 
the Hubbard U for u = 0. In this case A q = 0. 
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Figure 4.4: The gap parameter A c (in units of t) as function of t/U: the SCHA solution 

A SCHA _ e -2ntln(t/U)/U ^ &nd the Bethe angatz so l u ti n A^ a ' ~ e ~^t/U ^ both yalid j n 

the limit IT/* < 1 . 
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u/t=\ 
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u 

Figure 4.5: Charge gap A c and spin gap A s (in units of t) of the dimerized Hubbard 
model for U/t = 1 obtained within SCHA. We have used the value g s = 1 and the 
expressions g c = 1/ yl + U/ (irvp), f c ,s = ^Fyl (th^f)- 

For small f/ these results can be verified analytically. In order to solve the gap 
equations analytically we consider the case U > where g s = 1 and g c < 1, and 
restrict ourselves to the limit of small dimerization u. Using the analytical results for 
{4>l)tr and (0s)tr, given by Eqs. (J4.32J) - (|4.33J) . one obtains that for u — > the spin gap 
vanishes while the charge gap approaches a constant according to 

A c (m)-A c (0) oc m 4/3 , (4.40) 
A s (u) oc u 2/ \ (4.41) 

with cutoff-dependent prefactors (for details see appendix |FJ) . The exponent 2/3 that 
characterizes the opening of the spin gap is in accordance with the corresponding 
exponent of the dimerized antiferromagnetic Heisenberg chain up to a logarithmic 
correction in the prefactor [SB] . Since the Heisenberg model corresponds to the U — > oc 
limit of the half-filled Hubbard model, this indicates that the SCHA result (14.41)1 is 
exact and persists even in the strong-coupling regime U/t^$> 1. For u > u* the behavior 
of the gaps is changed to jHlj 

A c (u) w A s (u) oc u 2 '^~ 9c \ (4.42) 

where the crossover value u* is defined by A s {u*) = A c (0). The u — U phase diagram 
is shown in Fig. 14.61 In Fig. 14.71 we show A c (u) — A c (0) and A s (u) as function of u for 
U/t = 1 as obtained from the numerical solution of the gap equations. For comparison, 
the analytical results of Eqs. (|4.4())1 and (|4.41jl are also shown. 
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Figure 4.6: Phase diagram of the dimerized Hubbard model, u* 
from a dimerization gap to a correlation gap. 
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Figure 4.7: Charge gap A c (u) — A c (0) and spin gap A s (u) (in units of t) of the 
dimerized Hubbard model for U/t = 1 obtained within the SCHA. The straight lines 
are the analytic results of Eqs. 1)4.40)1 and 1)4.41)1 valid for u < u* . Note that A c (0) is 
subtracted from the charge gap in order to highlight the power law behavior. 
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4.2.2 Finite systems 

As in the case of spinless fermions, for a finite system it is not possible to simply 
replace the Klein factors by their eigenvalues since the terms proportional to the spin 
and charge current, J c and J s , in the Luttinger Hamiltonian (|4.4J) do not commute with 
the F's. One may, however, decouple the Klein factors from the bosonic fields using a 
variational ansatz. We introduce the Klein Hamiltonian 

h b cs b c b s _ iB cs L(F^F L ^ + F^Fli) + BcLF^F^F^Fli 
+B S LF+^ F+ F Rl F L i + h.c. 

+^{v c gcJc+v s gsJs), (4-43) 

where B c , B s and B cs are now variational parameters to be determined self-consistently. 
jjBcsBcBs j g Q £ form of a tight-binding Hamiltonian for a particle moving on a 2d 
lattice in a harmonic potential. The explicit representation of this Hamiltonian in the 
\J C , J s ) is given in appendix |Ej The minimum condition for the variational energy is 
equivalent to the following substitution in the backward and Umklapp terms 



F+F+F 7 F p e ±ifca ^ -> (F+F+F 7 F p )e ±i2 ^ Q + F+F+F 7 F p (e ±l2 ^ a ) (4.44) 
and to the substitution 

_p+_p^ e ±«v / 2</'aT«v / 206 _y ^_p+ jr^ e ±iV24>aTiV24>b _|_ JT+ jp^ (^ e ±iV2<f>a=FiV2<j>b^ (4.45) 

in the Peierls term, where a,/3, 7, p G {R ],L |, Jf2 i,L |} and a,b & {c, s}. As a 
result, instead of replacing the products of Klein factors by their eigenvalues as in the 
thermodynamic limit we now have to replace them by their expectation values with 
respect to the ground state of }j^, csBcBs divided by the system size. The details of 
the calculations are given in appendix EI The Klein Hamiltonian ()4.43|) can only be 
diagonalized numerically. In the bosonic sector one ends up with a sine-Gordon type 
model like in the thermodynamic limit. In order to evaluate the mean field parameter, 
in the framework of the SCHA, the quantity e$ which enters the gap equations (and 
which is explicitly given in Table l4~Tj) is replaced by e° r , the ground state energy of 
jjB cs B c B a (-[jyjjg^ by £ ( see appendix EJ. Apart from this modification Eqs. (|4.28j) - 
(|4.29j) remain unchanged. 

We consider a system with both spin and charge gap, i.e. ji^O. For large systems 
the parameters B cs , B c and B s become size-independent, with the consequence that the 
kinetic energy in H BcsBcBs dominates the confining potential. While approaching the 
thermodynamic limit, the e° r , is given by the minimum of e(k^, fcj), see Eq. (|4.24|) . and 
the expectation values of the Klein factors converge to the corresponding eigenvalues 
as can be seen in Fig. | 
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Figure 4.8: The derivative of e° r with respect to the variational parameters B c , B s , B cs 
needed to find the solution of equation (J4.28J) . ()4.29|) . for u ^ 0. With increasing L 
they converge to the values corresponding to the thermodynamic limit given in Table 
O de°JdB c = 2, de°JdB s = -2, de° tr /dB cs = -4. 
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Figure 4.9: The gap parameter A c (in units of t) as a function of the dimerization u 
for different system sizes compared with the thermodynamic limit. 



4.3 Ionic Hubbard model 
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Figure 4.10: The gap parameter A s (in units of t) as a function of the dimerization 
u for different systems sizes compared with the thermodynamic limit; A s has to be 
compared to the finite size energy Vp/L, which is e.g. Vp/L — 2 • 10~ 3 t for L = 1000. 



In Figs. 14.91 and 14.101 we show the dependence of the gap parameters A c and A s 
on the system size. Finite size effects are most pronounced in the small u region that 
corresponds to the Mott insulating phase. In particular the spin gap parameter A s 
approaches a finite value for all finite systems for u — > 0, while it goes to zero in the 
thermodynamic limit. 



4.3 Ionic Hubbard model 

The Hubbard model with an additional alternating on-site energy modulation has been 
used to study the organic mixed-stack charge transfer crystals with alternating donor 
and acceptor molecules [HHl EHl E3 EH] or to understand the ferroelectric transition in 
perovskite materials such as BaTiC>3 [BH] or KNb03[7U]. The model is described by 
the Hamiltonian 

H =-t E« c ^ + h - c -) + u E n n n n + A E(- 1 )' c i c ^> ( 4 - 46 ) 
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where A is the amplitude of the staggered potential. At U — and for A > 0, the 
model corresponds to a conventional band insulator with a band gap 2A. The one- 
dimensional half-filled Hubbard model without an alternating potential (A = 0) and 
with U > describes a correlated insulator with vanishing spin gap. 

In the atomic limit, t = 0, and for U < 2A, every second site of the lattice with 
on-site energy — A is occupied by two electrons while the sites with energy A are empty. 
The energy difference between the ground state and the highly degenerate first excited 
state is 2 A — U. For U > 2 A, each site is occupied by one electron and the energy 
gap is U — 2A. Within this simplified picture one obtains a critical point U C (A) = 2A, 
for which the excitation gap disappears. These considerations suggest that the system 
will be in two quantitatively different phases in the limits U <^ A and U 3> A, with a 
quantum phase transition from a band insulator to a correlated insulator in between. 

The early numerical [ZU E2j and analytical HI] results reported that at T = 
and for fixed A a single phase transition occurs if U is varied. Using bosonization 
Fabrizio et al. [73] predicted a two — transition scenario for the ground state phase 
diagram. For U < U& the system is a band insulator with finite charge and spin gaps. 
The first critical point U c \ is an Ising critical point, where the charge gap vanishes. The 
intermediate phase, for U c \ < U < U C 2, is a spontaneously dimerized insulator phase, 
in which the bosonic spin and charge gaps are finite. This phase is characterized by a 
non-zero expectation value of the dimerization operator 



The second transition, at U — U C 2, is of the Kosterlitz-Thouless type, and the system 
goes over into a correlated insulator phase with a finite charge gap and a vanishing spin 
gap. Various attempts based on numerical tools have been made to verify this scenario, 
without leading to a definite conclusion. In particular, exact diagonalization [Tol I77j . 
valence bond techniques [7H], quantum Monte Carlo [ZH] and the DMRG |%U ] IHT ] IH2] 
were used with different results in favor of one or two critical points. 
In momentum space the Hamiltonian (|4.46J) reads 



where iconic is the contribution due to the alternating on-site energy modulation 




(4.47) 






(4.49) 



fe<7 



which is referred to as the "ionic" term. In bosonized form it reads 
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+ A [ L dxF+Fue-^-^e-^*-^ + h.c 
At half filling, N c = N s = 0, the ionic contribution is 

H ionic = A Tdx {F+ T F LT e-^ + ^) + F^F Ll e~ iV ~ 2 ^-^ + h.c.}, (4.50) 
Jo 

where A = A/27ra. Similar to the case of the Peierls-Hubbard model, we express the 
Klein factors in terms of the operators = F^F^ and Ai = F^F^. We choose 
a basis in which these operators are diagonal, and using Eqs. (|4.1Hj) and (|4.1fij) . we 
replace the Klein factors in Hi, H 2 and H ionic by their eigenvalues to obtain 



Hi = U dx { e ^-h) e -^V2<t>s + h c _} ; 
Jo 

rL 

H 2 = U dx { e *(fcT+fc|) e -*2v / 20 c + hx _| ; 
Jo 

H ionic = A / dx { e ^r e -iV2{<t> c +^) + e ik le -iV2(^-^) + h c y 
Jo 



(4.51) 



(4.52) 



(4.53) 



As in the case of the dimerized Hubbard model, the Hamiltonian of the ionic Hubbard 
model separates into different sectors of purely bosonic Hamiltonians, labeled by fc-f 
and ki- Using Majorana fermions ^3] to represent the Klein factors, only one sector 
can be obtained and again the continuity is lost. 

4.3.1 Self-consistent harmonic approximation 

Since the nonlinear terms of the ionic Hubbard model cannot be treated exactly we 
again employ the SCHA as in the case of Peierls Hubbard model. In the thermody- 
namic limit we replace the Klein factors by their eigenvalues and introduce the trial 
Hamiltonian 

= E r ? \-( 9 M 2 + v a9a (d x e a ) 2 + -^A , (4.54) 

^ s Jo 2tt [g a v a g a J 

which provides us with a variational estimate for the ground state energy 

E = (H ) tI + (Hi) tT + (H 2 ) tr +(H lonic ) tr , (4.55) 
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B c > B s , range 
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1 
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( TD TD TD \ 


< B cs < 2^B C B S 





IT 


— 2B C — 2B S 


2\JB C B S < B cs < 2B C 


arccos + n 


arccos Hr — ix 


-2B C + 2B s -§^ 


2B C < B cs 


TX 


IT 


2B C + 2B S — AB CS 



Table 4.2: Minimum eo of e(k^ , /q) = 2B C cos(k^ + /q) + 25, cos(k^ — /q) + 2i? cs (cos /c-f + 
cosfcj) which is used in the gap equations ()4.6H) and ()4.62|) . The first line is only 
relevant for A = (Hubbard model without on-site energy modulation), whereas the 
second and the third lines are relevant for A > 0. 



where E tr is the ground state energy of H tv , and 

e(/cj, ki) = 2B C cos(/c-f + /q) + 2B S cos(/c-f — + 2£> cs (cos k^ + cos k\), (4.57) 

with 

B c = f/ e - 4( ^ >tr , (4.58) 

B s = f/e" 4 ^, (4.59) 

B cs = Ae- (fi "e- ( *" . (4.60) 

Minimizing the variational ground state energy with respect to A c and A s yields the 
gap equations 

4B ^wb~ - ticsw^, (4.61) 



2nv c g c dB c dB, 

^- = -v.* -3°%-, (4.62) 
2-Kv s g s dB s dBcs 

where eo is the minimum of e(k^ki) with respect to k^ and k± (see Table l4~2*j) . De- 
pending on the values of B c , B s and B cs , different physical situations arise: 

1) The solution eo = —2B C — 2B S corresponds to the Hubbard model without on- 
site energy modulation, A = 0. The system is a Mott insulator with a vanishing 
spin gap, as in the Peierls-Hubbard model with u = 0, see Fig. 14.31 and Eqs. 
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2) The solution eo = 2B C + 2B S — 4B CS corresponds to a band insulator, where the 
charge and spin gap decrease with increasing U . In this case the gap equations 
l4~oTl) - (l4~7)2l read 



A 2 

c 



2irv c g c 
A 2 

s 

2nv s g s 



AB CS - 8B C , (4.63) 
AB CS - 8B S . (4.64) 



As shown in appendix |Fl this case corresponds to the limit A/U ^> 1 , and one 
obtains the following dependence of the gap parameters A c and A s on the ionic 
potential A: 

A c oc A 2/(3 " 9c) , (4.65) 



A s oc A 2/(3 ~ 9c) . (4.66) 

This is the same result as obtained for the Peierls-Hubbard model with u replaced 
by A. For U = 0, i.e. g c = 1, Eqs. ()4.65|) - (j4.66J) yield A Cj s oc A in accordance 
with the known solution in the non-interacting limit. 

3) The solution e = —2B C + 2B S — B 2 S / B c corresponds to the spontaneously dimer- 
ized insulator phase where the dimerization operator has the expectation value 



B- 2 



(D) ~ sin(/c T ) + sin(^) = -2W 1 - ^ j . (4.67) 
In this case the gap equations are 

A 2 r>2 

8fi c -2-f, (4.68) 



2wv c g c B c 

A 2 n2 

= -8Bs + 2-^, (4.69) 
2nv s g s B c 

which are solved in appendix |FJ The result is 

A C (A) - A c (0) oc A 4 , (4.70) 



A a (A) oc A 2 , 



(4.71) 
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Figure 4.11: The mean field parameters A c and A s (in units of t) as function of A, 
the staggered potential, in the thermodynamic limit. At A ~ 5 • 10~ 4 t we find a 
discontinous jump of the parameters. 



which is similar to the behavior of the Peierls-Hubbard model, Eqs. (I4.40J) - ()4.41J) . how- 
ever with different exponents, 4 instead of 4/3 for the charge gap, and 2 instead of 
2/3 for the spin gap. This different behavior can be understood in the strong coupling 
limit U — > oo where in second order perturbation theory the Peierls term yields an 
alternating contribution to the Heisenberg exchange coupling J = 4t 2 /U, while for 
the ionic Hubbard model the strong coupling limit is the antiferromagnetic Heisenberg 
model with uniform coupling. In this case the opening of the spin gap is associated 
with higher than second order processes in the strong coupling limit. 

Fig. 14.111 shows the dependence of A c and A s on A for fixed Hubbard interaction 
U/t = 1. Contrary to the case of the Peierls-Hubbard model there is no smooth 
crossover between the small and large A region but a discontinuous transition at some 
intermediate value A*. In Fig. 14.121 the ionic parameter A = 10~ 3 is kept fixed and A c 
and A s are plotted as functions of U. For small values of U spin and charge gaps are 
almost equal and decrease with increasing U, as highlighted in Fig. 14.131 At U/t ~ 1.05 
there is a discontinuous transition to a phase where A s goes to zero while A c increases 
strongly, as can be seen in Fig. 14.141 

Connecting the values of A where the discontinuous transition occurs for different 
values of U one obtains the phase boundary between a Mott insulator (A s — > 0) and a 
band insulator (A c ~ A s ) as shown in Fig. 14.151 We find only one phase transition, and 
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Figure 4.13: The mean field parameters A c and A s (in units of t), as function of U, in 
the thermodynamic limit. Note the different scales compared to the previous figure. 
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Figure 4.14: The mean field parameters A c and A s (in units of t), as function of U in 
the thermodynamic limit. Note the different scales compared to the previous figures. 
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Figure 4.15: Phase diagram of the ionic Hubbard model. The dots correspond to the 
discontinuity of A c and A s , as for example shown in Fig. 14. Ill for U/t = 1. 
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not two, in contrast to the phase diagram proposed by Fabrizion et al. [75]. However, 
in our bosonization approach, using SCHA to estimate the charge and the spin gap, 
we are limited to small values of U and A. 
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4.3.2 Finite systems 

To calculate the gaps of the ionic Hubbard model for a finite system we apply the same 
procedure as in the case of the Peierls-Hubbard model. We will find, unexpectedly, 
strong finite size effects near the phase transition. At present it is still not clear to us 
whether this behavior is an artefact of the SCHA or if it has physical reality. 

We decouple the Klein factors from the bosonic fields introducing the Klein Hamil- 
tonian 



where B CJ B s and B cs are variational parameters. H tI cs c s is again of the form of a 
tight-binding Hamiltonian for a particle moving on a 2d lattice in a harmonic potential. 
The explicit representation of this hopping Hamiltonian in the \J C , J s ) basis is given in 
appendix |E| for the Peierls-Hubbard model. For the ionic Hubbard model one has to 
replace ±iB cs by B cs to obtain the corresponding matrix elements. As a result of the 
variational procedure, the Klein factors are replaced by their expectation values with 
respect to the ground state of the Klein Hamiltonian. Similar to the case of the Peierls- 
Hubbard model, in the framework of the SCHA, the quantity eo which enters the gap 
equations (and which is explicitly given in Table is replaced by e° r , the ground state 
energy of f^ csBcBs divided by L, which can be calculated only numerically. Fig. 14.161 
shows the derivatives of e° r with respect to B c , B s and B cs as function of the system 
size for fixed model parameters U/t = 1 and A/t = 0.1. For small L the derivatives are 
close to zero while with increasing L they converge to the values corresponding to the 
thermodynamic limit given in Table l4~2l de® T /dB c = 2, de® T /dB s = 2, de® T /dB cs = —4. 
The scale 1 / \fL is chosen in order to enlarge the crossover region. 

In Fig. 14.171 we show the gap parameter A c as function of A and fixed Hubbard 
interaction U/t = 1 for different system sizes. While for large values of A, correspond- 
ing to the band insulator phase, there are nearly no finite size effects, for small A, 
corresponding to a Mott insulator phase, even for L = 5000 the values have not yet 
converged to the thermodynamic limit. The transition between the Mott insulator and 
the band insulator phase appears continuous for small values of the system size L but 
becomes sharper with increasing L and finally turns over into a discontinuous one. The 
same observation holds for the spin gap parameter A s , which is displayed in Fig. 14.181 
for the same set of parameters as before. In particular for A — > the gap parameter 
A s approaches a constant for finite systems while in the thermodynamic limit it goes 
to zero according to A s oc A 2 . 

In order to show in more detail how the discontinuous transition emerges from the 
continuous one we plot the gap parameters A c and A s for several values of L in the 



r B cs B c B s 
tr 



B cs L(F^F L i + Fj^Fn) + B c LF^F^F Ll F L ^ 
+B S LF+ F+F Rl F L i + h.c. 

+ ^y(v c g c J 2 c +v s g s J 2 s ), 



(4.72) 
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Figure 4.16: The derivative of e° r with respect to the variational parameters B c , B s , B, 
needed in the equations (|4.61|) and 1)4. 62j) . 
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Figure 4.17: The charge gap parameter A c (in units of t) as a function of A for different 
values of the system size L compared with the thermodynamic limit. 



4.3 Ionic Hubbard model 



59 



0.1 C 



< 

< 



le-05 



U/t = 1.0 

1 1 1 — I I I I I | 



T 1 1 — I I I I IJ 



0.01 r 



" Q ...0...fl-E-E-Q-"; 

o.ooi r_ t n^r 



A 



0.0001 \r * 



L = 


2000 


■-B-- 


L = 


5000 




L = 


00 





A' 



4' 



J I I 



1 



J I I 



0.001 

A/t 



0.01 



Figure 4.18: The spin parameter A s (in units of t) as a function of A for different 
values of the system size L compared with the thermodynamic limit. 
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Figure 4.19: The charge parameter A c (in units of t) as a function of A for different 
values of the size system size L in the region where the transition from a band insulator 
to a Mott insulator appears. 



region where the transition occurs in Figs. I4~TT?1 and I4".20l respectively. For U/t = 1, 
the discontinuity appears somewhere between L = 5000 and L = 5500 and becomes 
more pronounced while L is increased further. 



In Fi gs . 14 . 211 and I4~2*2l we show A c and A s , respectively, as function of the Hubbard 
interaction U/t for fixed value of the ionic parameter A/t = 10~ 3 and different system 
sizes L. For this small value of A/t there are still quite pronounced finite size effects for 
L = 1000. This indicates that numerical methods which are devised for finite systems 
like exact diagonalization or DMRG would probably fail in the small A regime. 




Figure 4.20: The spin parameter A s (in units of t) as a function of A for different 
system size L in the region where the transition from a band insulator to a Mott 
insulator appears. 
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Figure 4.21: The charge parameter A c (in units of t) as a function of the Hubbard 
interaction U/t for different system sizes compared with the thermodynamic limit. 
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Figure 4.22: The spin parameter A s (in units of t) as a function of the Hubbard 
interaction U/t for different system sizes compared with the thermodynamic limit 
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4.4 Conclusions 

In this chapter we considered two different one-dimensional Hubbard models, the 
Peierls-Hubbard and the ionic Hubbard model. In the bosonized version of these mod- 
els two kinds of bosonic fields have to be introduced which are associated with the 
spin and the charge degrees of freedom, respectively. The backward and the Umklapp 
scattering yield nonlinear terms in the spin and the charge sector, respectively, while 
the Peierls distorsion and the ionic potential introduce a coupling between the spin 
and the charge sectors. This nonlinear terms lead to the opening of spin or charge gap, 
depending on the model parameters. 

Our first goal was to stress the importance of Klein factors in bosonized Hamiltoni- 
ans with nonlinear perturbations. In the study of such types of Hamiltonian, one has 
to keep in mind that the Klein factors which accompany the nonlinear terms do not 
commute with the total spin and charge currents, J CjS , and therefore Klein factors and 
current operators cannot acquire a fixed value at the same time. In the thermodynamic 
limit the ground state of the gapped system is a superposition of many states with dif- 
ferent J c and J s . In this situation it is possible to choose a fixed phase for the Klein 
operators. The bosonized Hamiltonian is then the conventional sine-Gordon Hamil- 
tonian. For finite size systems the Klein factors can be replaced by their expectation 
values with respect to the ground state of a trial Hamiltonian, and the resulting bosonic 
Hamiltonian is also of a sine-Gordon type, however with size-dependent couplings. 

Using a self consistent variational scheme, we decoupled the bosonic Hamiltonian 
from the Klein Hamiltonian and we calculated the gap parameters A c and A s as a 
function of the strength of the perturbation (the dimerization u in the case of the 
Peierls-Hubbard model, and the ionic potential A in the case of the ionic Hubbard 
model). In both cases we found a transition from a Mott insulator phase with vanishing 
or very small spin gap to a band insulator phase where spin and charge gaps are of 
the same order of magnitude. The main difference between the two models is that 
the transition between the Mott and the Peierls phase is continuous in the case of 
Peierls-Hubbard model and discontinuous for the ionic Hubbard model. 



Chapter 5 



Screening in low-dimensional 
electron systems 

In order to extend the model calculations of the preceding sections towards a more 
realistic description of quasi one-dimensional systems we consider in the following a 
one-dimensional interacting system, characterized by a Luttinger model, coupled to a 
three-dimensional environment, and calculate the spin and the charge susceptibility. 
The motivation is an old discussion on the role of electronic correlations for the metal- 
insulator transition in VO2. 

5.1 Phase transition in VO2 

A detailed analysis of electronic structure of vanadium dioxide (V0 2 ) is given in [SHI Ell- 
in the following we recall some of the theoretical and experimental findings. Vanadium 
dioxide in the paramagnetic metallic phase above T c = 340 K has a rutile crystal 
structure, each vanadium atom being located at the center of an oxygen octahedron. 
The fivefold-degenerate d levels of the V 4+ (3<i 1 ) ion are split into doubly degenerate e g 
levels and triply degenerate ti g levels in the octahedral crystal field. The e g orbitals 
are strongly hybridized with the O 2pa orbitals and have a large bandwidth, and are 
pushed far away from the Fermi surface. The t2 9 levels are split into the du and n* 
levels by the orthorhombic component of the tetragonal crystal field. Thus the du and 
7r* bands are situated at the lowest energies around the Fermi level. According to LDA 
calculations |SSJ E3j the d\\ bands show a strongly one-dimensional dispersion, whereas 
the ir* bands are three-dimensional. These bands play the crucial role in electronic 
(transport and magnetic) properties VO2, and will be further studied in section I5~^l 

In the insulating non-magnetic phase below T c , VO2 is distorted to a monoclinic 
crystal structure, involving a pairing of vanadium ions V 4+ along the c r axis. Because 
of the change in the V-0 hybridization, the energy of the hybridized ir* band rises 
above the Fermi level and it becomes empty. Furthermore, the du band is split into 
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two bands. A schematic energy diagram of the 3d bands around the Fermi level for 
VO2 is shown in Fig. 15.11 

There is an old discussion concerning the driving mechanism for the metal-insulator 
transition in this material and the role of electronic correlations. Some authors argued 
that electronic correlation effects are in VO2 less pronounced than in the other vana- 
dium oxides, like V2O3. Goodenough pointed out that the metal-insulator transition 
of VO2 may be explained by considering the change in crystal structure, which leads to 
the splitting of the d\\ band and the rising of the it* bands [H3]. Modern LDA jEHlEI] 
and molecular dynamics [EH] calculations support this point of view, although the split- 
ting of the d\\ bands is underestimated, so that they predicted a metallic instead of an 
insulating state. 

Further experimentals, however, demonstrated that electronic correlations play an 
important role also in VO2. Pouget et al. [S3 EE! performed nuclear magnetic resonance 
(NMR) measurements under uniaxial stress in VO2 and Vi^Cr^C^. They found for 
Vi_ x Cr x 02 three different insulating phases (Mi, M2 and T). The Mi phase is identical 
to the insulating phase of pure VO2 at normal pressure, while in the M2 and T phases 
magnetic moments were observed. The appearance of the magnetic moments was 
considered as a clear sign for correlation effect in the insulating phase. 

Zylbersztejn and Mott emphasized the important role of the it* electrons in the 
metal-insulator transition. They argued that, in the insulating phase, the hybridization 
between the the tt* and O 2pa orbitals increases due to the crystal distortion so that 
the it* band rises in energy and becomes empty. Therefore the correlation energy of 
the d\\ electrons becomes enhanced leading to a splitting into the lower and the upper 
Hubbard bands. They argued that the metallic phase is adequately described within 
the framework of band theory, since the interaction of the d\\ electrons is screened by 
the 7r* electrons. They estimated the effective intra atomic Coulomb repulsion in the 
metallic phase from the Stoner enhancement of the magnetic susceptibility as 

$=inbr M (5 - 1) 

yielding £/" eff = 0.13 eV. Here xf is the Pauli spin susceptibility and Md\\ = 5.25 
states/eV is the density of states per spin at the Fermi surface in the d\\ band. 

Shin et al. [EH! measured the vacuum-ultraviolet reflectance and photoemission spec- 
tra of VO2 in order to investigate the 3d bands and electron-correlation effects. The 
splitting energy between the two d\\ bands in the insulating phase was found to be about 
2.5 eV, bigger than the value obtained from the cluster band calculation (§0], namely 
0.5 eV. They argued that in the insulating phase the band splitting of the it* — dn and 
dn — d\\ bands cannot be explained by the crystal distortion alone, correlation effects 
playing an important role also. They estimated the correlation energy in the insulating 
phase of the d\\ band as about U(d\\,d\\) ~ 2.1 eV. From the photoemission data they 
concluded that the correlation energy of the dn electrons must not be neglected even 
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Figure 5.1: Schematic energy diagram of the 3d bands around the Fermi level for VO2 
|89j . The energy of 0.2 eV is obtained by the UPS spectra. The energies of 0.5 and 1.1 
eV are obtained by ultraviolet reflectance spectra, while 0.65 eV is a value from the 
infrared absorption spectra jHIj 



in the metallic phase, and U(d\\,d\\) was estimated to be not less than 1.3 eV. 



5.2 The model 

Following the ideas of Mott and Zylbersztejn that in the metallic state of VO2 the 
electron-correlation effects in the one- dimensional d\\ bands are screened by the three- 
dimensional 7i* electrons, we consider a one-dimensional, interacting Fermi system, 
which is coupled to a three-dimensional environment through a local interaction. The 
aim is to see how the presence of the three-dimensional system affects the one-dimensional 
system. We calculate the spin and charge susceptibility, with the result that the spin 
susceptibility is not changed due to the coupling to the environment, whereas the 
charge susceptibility shows screening. 
We consider the Hamiltonian 

H = -f^Luttingcr + -^Ott* + Hd-n* , (5-2) 

where iiLuttinger is the Hamiltonian corresponding to the Luttinger model, 

= J2 e ^<>L (5-3) 

corresponds to a band of non-interacting it* electrons with a three-dimensional disper- 
sion e n (k), and 

H d-** = d L d k-wKt'K+w' ( 5 - 4 ) 

kpq,CTcr' 
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describes coupling between tt* and d\\ electrons via an interaction of strength V; N is 
the number of lattice points (in the three dimensional lattice). In order to make a con- 
nection between our calculation and the discussion on VO2, we use the notation d + , d 
for the creation and annihilation operators of d\\ (one-dimensional) electrons, and the 
notation tt* + ,tt* for the creation and annihilation operators of tt* (three-dimensional) 
electrons. A similar Hamiltonian has been studied by Paquet and Leroux-Hugon [§2*] . 
In this work, electronic correlations in the one-dimensional du band were taken into ac- 
count in terms of a spatially inhomogenous mean field theory. In contrast, we describe 
the du electrons as Luttinger liquid. 



5.3 Charge and spin susceptibility 

The dynamic susceptibility \ is a measure of the system's response to an external field. 
Generally, the dynamic susceptibility is given by a retarded commutator of the form 

XaaM =z0(t)([A(q,t),A(-q,O)]). (5.5) 



For the spin susceptibility A is the spin density operator, er(q) = P|(q) — Pj,(q), while 
for the charge susceptibility A is the charge density operator, p(q) = p-f(q) + Pj.(q). 
The density operators p^j are given by Eq. (|2.8|) where for the general case the one 
dimensional momenta have to be replaced by the three dimensional vectors. 



5.3.1 Luttinger model 

For the calculation of the charge susceptibility in the Luttinger model we express the 
charge density in terms of left and right charge densities, 

p(q) = Pr(q) + pM- 

Corresponding to this separation in left and right movers the correlator Xpp can De 
expressed as 

Xpp = iO(t)(\p(<l, t),p(-q, 0)]) = Xll + Xlr + Xrl + Xrr- (5-6) 

By taking the time derivative one finds for all four components Xvy an equation of 
motion of the type 

-idtXm = S(t)([p v (q,0), Pj (-q,0)}) - ie(t)([\p v (q,t), H], p y (-q,0)]). (5.7) 
For the "g-ology" Hamiltonian the commutators [p v , H] are 



[p L , H] = qv F (l + 74|| + 74±)pL + <Ff(-7i|| + 72|| + 72±)p«, 



(5.8) 
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[p R , H] = -qv F (l + 7 4 || + 74±)pr - <?M-7i|| + 72|| + l2±)pL, (5.9) 

where ji = gi/(2%VF)- With this simple form of the commutators one can calculate 
the susceptibility explicitly, and the final result is 



t 9cVF -{qv, 



XpP (q,uj) = 2NaN ld {e F y c ~* ^ = x %, (5.10) 

v c io l - [qv c y pp 

where Mw^f) = (t^f) -1 is the density of states per spin at the Fermi surface, a is 
the lattice constant, and v c and g c are the Luttinger parameters in the charge sector 
as given in Eqs. (|2.52|) and (|2.54j) . 

In a similar way, the spin density can be separated into left and right spin densities 

<7(5) = 0*(?) + ox(s)- (5-11) 

Following the same steps as for the charge susceptibility one obtains for the spin sus- 
ceptibility 

X **(q,u) = 2NoM ld (e F ) 9sVF ~^ s) \ 2 , (5.12) 

where v s and g s are the Luttinger parameters in the spin sector as given in Eqs. (|2.51j) 
and fl233). 

5.3.2 Luttinger model coupled to 3d electrons 

Now we calculate the charge susceptibility of a one-dimensional electronic system which 
is coupled to a three-dimensional environment. The model is described by the Hamil- 
tonian ()5.2j) 

H = -^Luttinger + H 07T * + Hd- n * ■ (5.13) 

We consider in the following only interactions with small momentum transfer; in this 
case, in terms of density operators, the Hamilton H d - n * reads 

q 

where q is the momentum transfer with q = (gy, qj_). Since the total charge (and 
spin) density is a sum of three contibutions, p = Pl + Pr + Pn*, then the correlator \pp 
is 



X PP = XLL + XRL + XLR + XRR + XLtt + XRn + Xtttt + XttL + XttR- 



(5.15) 
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To calculate the correlators Xrryi we use again Eq. (|5.7j) . this time for 77, 7 e {L, i?, 7r*}. 
In addition we need the commutators 

([^(q,0), P7r .(-q,0)]) = ^(K* CT (k+q))-K* CT (k))) « -^(e^Vv^-q, (5.16) 

ko- 



([[p ff .(q, t), ffftr*], Pr,(-q, 0)]) « -v| d • q([^.(q, t), p„(-q, 0)]), (5.17) 



([b^(q,*),-^d-7r*],^(-q,0)]) = 

-2tf 3d (e F )v 3 F d ■ q^V([ PL (q, t) + p R (q, t), p v (-q, 0)]), (5.18) 



[p L (q,t),H d ^] = -[p R (d,t),H d ^] = ^^Mp^(q,t). (5.19) 

where N± is the number of chains , M-sd^F) is the density of states per spin and per 
volume at the Fermi surface, and V is the volume. After some calculations we obtain 
the following expressions for the charge and spin susceptibility: 

v _ X%-2X%VX«-/N + Xv _ 



X<r* = X% + X«; (5-21) 

where x]f P (xlt) is the charge (spin) susceptibility of the Luttinger model given by Eq. 
1)5.1 Oj) (Eq. ()5.12j) ) with Na = LN±, and Xn* is the susceptibility corresponding to 
the free tt* electrons. The results ()5.20|) and ()5.21|) are equivalent to the RPA results 
obtained for a similar model with an electron-electron interaction characterized by a 
small momentum transfer, without applying the bosonization procedure. 

5.3.3 Dispersion of the charge excitations 

The charge excitation spectrum is given by equation 
1 



x PP (q, 



0, (5.22) 



which implicitly defines the dispersion relation uj = uj(q). In the following we analyze 
the dispersion of the charge excitations: We want to see how the linearized dispersion of 
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a clean one- dimensional system is modified due to the presence of the three-dimensional 
electrons. In this case Eq. (15.22)1 is leads to 



q\\)V X Au, q)V/^ = 0, (5.23) 



where Xopi^i Q\\) 1S the dynamical susceptibility of the Luttinger model model, compare 
Eq. ((03), 



X >,?||) = ZNaJsrjm 2 { T e) \ 3 , (5-24) 

and x'LO^q) is the dynamical susceptibility of the free three-dimensional fermionic 

system, 

X^(w,q) = 2AT 3d v/ V ^ M q \ , (5.25) 

where A/m, v F d are the density of states per spin and per volume and the velocity of the 
7r* electrons at the Fermi surface and (...)fs is a Fermi surface average. We separate the 
three-dimensional momentum q in two components, one parallel to the one- dimensional 
chain, q\\, and the other perpendicular to it, q_i_. If we insert Eqs. (15.24)1 and (J5.25)) 
into Eq. (|5.23j) . the dispersion is found as the solution of the following equation: 



r 3d 

.2 /„ „. \2 



w - (q\\v c ) 



d N v r \u + v^-q/ FS 



0. (5.26) 



In the limit u Cv^ -q the three dimensional charge susceptibility can be approximated 
by its static value 

Xir* {u, q) ~ Xtt* (0, q) = 2Af 3d V 
and the dispersion of the one- dimensional charge excitations is given by 



u = d^ |M /l-4V>M*A^f^— ■ (5-27) 

JX v c 

In this case the dispersion maintains its one-dimensional linear characteristic. The only 
effect of the coupling to the three-dimensional environment is a renormalization of the 
Luttinger parameters 



4V>M^ M ^^, (5-28) 

JX V r 



J C 
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which is equivalent to a decrease of the interaction strength in the charge channel. 
In the opposite limit, to ^> ■ q we can expand in powers of • q/u>, and obtain 



• q 

,3d 



i«) 2 (g +qi) 



w -j-vf ■ q/ FS 3 uj 2 
In this situation the dispersion is given by the equation 



■ aV ^g c v F ( j, 
N 3v c 



u A - uj\q\\v L ) 2 + {q\\v c ) 2 V 2 N ld Nzd 
with the following solutions: 



v 



r) 2 «+qi; 



/ ld- / V3 ( i 



(?|| u c) : 



(5.30) 



(5.31) 



(5.32) 



ul = V 2 M ld M^-^{vf) 2 {q 2 + q 2 J. (5.33) 

In the weak coupling limit the solution u 2 , is not compatible with the assumption 
oo v^ d -q. To analyze the solution uj 2 ^ we must distinguish between two limiting cases. 
If the parallel component of the momentum is much larger than the perpendicular 
component q\\ 3> q±, then again the effect of the tree-dimensional electrons is only to 
renormalize the Luttinger parameters and the dispersion remains approximately, linear 



\ 



v c 



3v r N 



(5.34) 



A more interesting situation arises in the opposite limit, q\\ <C q±, when the charge 
dispersion looses its linearized characteristic and the one dimensional modes are mixed 
with the three-dimensional modes through the q± component, i.e. 



U>4 



(Q\\ v ) 



3v r N 



(5.35) 



with 



a 



vf ± ) 2 Ag c v F aV 
vl 3v r N 



V 2 Af ld Af 3d 



V 2 



, Id ,3a! 
t F t F 



where we have assumed that V is roughly one order of magnitude smaller than the 
Fermi energies. 



5.4 Discussion 
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Figure 5.2: The spin susceptibility of the one-dimensional Hubbard model. Within 
RPA one finds Xs(U) = x s (0)/(l — UAfid), where Mid — (^f)" 1 - The data points are 
numerical results jHS] obtained with exact diagonalization of a Hubbard chain with 6 
electrons on 9 lattice sites (n is the electron density). 

5.4 Discussion 

In order to see how the three-dimensional electrons affect the electrons in the one- 
dimensional d\\ band, we have calculated the charge and the spin susceptibility of the 
system; the coupling is introduced as a local Coulomb interaction of strength V. We 
found that the spin susceptibility of the dn electrons remains unchanged in the presence 
of the 7r* electrons. This seems to contradict |93j, where from the spin susceptibiliy 
in the metallic phase a smaller U e g ~ 0.13 eV has been estimated. The apparent con- 
tradiction is easily solved, however, when calculating the spin susceptibility in a Id 
Hubbard model numerically, and comparing with the RPA estimate. Fig. 15.41 cleary 
demonstrates that for a Stoner enhancement of about x/Xo ~ 3 the RPA drastically 
underestimates the interaction strenght which is necessary to obtain such an enhence- 
ment. 

Concerning the charge susceptibility, we found that it is changed in the presence 
of the 7r* electrons. To make the effect more visible, we analyzed the dispersion of the 
collective charge excitations. When the dynamics of the tt* electrons is fast compared 
to the du electrons, we found that the major effect of the tt* electrons is to renormalize 
the charge velocity. When the Fermi velocity of the dn and ir* electrons are comparable 
in size, the dynamics of the it* electrons is on comparable time scales. In this limit 
collective modes can mix, so that the dn electrons loose their one-dimensional character. 
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Chapter 6 

Summary and outlook 



In this work we have investigated the low energy properties of different one-dimensional 
fermionic lattice models using the bosonization technique. We have attached much 
importance to a proper consideration of the Klein factors which are neglected or inac- 
curately treated in the majority of the literature. 

In order to deal with the nonlinear terms, arising from the dimerization or similar 
scattering processes, we used the self-consistent harmonic approximation (SCHA), a 
variational method where the nonlinearities are replaced by a quadratic potential. We 
have constructed an extension of the SCHA in order to account for the existence of 
Klein factors in bosonized Hamiltonians. 

As a first application, we have investigated a model of spinless fermions with 
nearest-neighbor interaction and a modulated hopping. For the infinite system, both, 
the value of the Luttinger parameter g c = 2 where the transition from a gapless to a 
gapped phase takes places, and the exponent 1/(2 — g) that characterizes the open- 
ing of the gap for small dimerization, are correctly obtained within the SCHA. In the 
thermodynamic limit replacing the Klein operators by a c-number is justified. When 
considering a finite system a more careful treatment is necessary. In this case the Klein 
Hamiltonian can be mapped onto a Mathieu equation which can be solved analytically 
in certain limits. Within our approach it turns out that the effect of the Klein factors 
is negligible as long as the dimerization gap is larger than the finite-size gap, vp/L. 

In addition, the finite-size formalism allows calculation of the size- dependence of the 
Drude weight. The Drude weight reflects the sensitivity of the system with respect to a 
change of boundary conditions and is related to the properties of the current operator J 
in the bosonized version of the Hamiltonian. In a finite system with an energy gap A the 
Drude weight is expected to be nonzero but exponentially small, D ~ exp(— const -LA). 
In our extended version of the SCHA we confirm the exponential behavior. 

As an application to fermions with spin, we studied the role of Klein factors for the 
Peierls- Hubbard and the ionic Hubbard model. In the bosonized version of these models 
two kinds of bosonic fields have to be introduced which are associated with the spin and 
the charge degrees of freedom, respectively. Nonlinear terms arising from backward and 
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Umklapp scattering lead to the opening of spin or charge gaps depending on the model 
parameters. Using the SCHA we have calculated the corresponding gap parameters 
A c and A s as a function of the strength of the perturbation (the dimerization u in the 
case of the Peierls-Hubbard model, and the ionic potential A in the case of the ionic 
Hubbard model). In both cases we found a transition from a Mott insulator phase 
with vanishing or very small spin gap to a band insulator phase with almost equal spin 
and charge gaps. The opening of the spin gap with increasing perturbation strength 
occurs according to a power law with universal exponents, 2/3 in the case of the Peierls- 
Hubbard model, and 2 in the case of the ionic Hubbard model. The main difference 
between the two models, however, is that in the Peierls-Hubbard model the transition 
between the Mott and the Peierls phase is continuous, while for the ionic Hubbard 
model a first order transition is observed, with a jump in A c and A s . As to the ionic 
Hubbard model we could not confirm the scenario proposed by Fabrizio et al. [ZS] who 
argued in favor of two phase transitions. However, in our bosonization approach we 
are strictly limited to weak coupling with regard to the Hubbard interaction U and 
to the ionic potential A, while Fabrizio et al. [ZH] used a strong coupling approach in 
order to obtain the second transition. While there is some evidence for the two phase 
transitions scenario from DMRG calculations in the case of a strong ionic potential 
(A/t = 0.5 in Ref. [S2] and A/t = 20 in Ref. [HI]), it is not clear whether this scenario 
persists to small values of A, e.g. A/t = 0.1 that we have used in our calculations. 

Concerning the role of Klein factors one has to recall that they do not commute 
with the total spin and charge currents, J CjS , and therefore Klein factors and current 
operators cannot acquire a fixed value at the same time. In the thermodynamic limit 
the ground state of the gapped systems is a superposition of many states with different 
J c and J s . In this situation it is possible to choose a fixed phase for the Klein operators. 
The bosonized Hamiltonian is then the conventional sine-Gordon-like Hamiltonian. In 
the framework of our variational scheme we decoupled the bosonic Hamiltonian from 
the Klein Hamiltonian. For finite size systems the eigenvalues of the Klein factors 
are replaced by their expectation values with respect to the ground state of the trial 
Hamiltonian, and the resulting bosonic Hamiltonian is also of a sine-Gordon type. 
The Klein Hamiltonian is of the form of a tight-binding Hamiltonian for a particle 
moving on a 2d lattice in a harmonic potential. The strength of the potential is 
inversely proportional to the length of the systems, i.e. for small systems the potential 
contribution to the ground state energy is relevant. On the other hand, when the size 
of the system approaches the thermodynamic limit, the confining potential becomes 
irrelevant. 

In an attempt to describe more realistic systems we considered a model where one- 
dimensional electrons are coupled to a three-dimensional conduction band via a local 
interaction, having in mind an application to embedded Peierls systems like vanadium 
dioxide, VO2. In order to see how the one-dimensional electron system is modified in 
the presence of the three-dimensional environment, we calculated the charge and the 
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spin susceptibility, and analyzed the collective excitations within RPA. We found that 
the spin excitations remain unaffected. In the charge channel, the collective modes 
of the Id and 3d subsystems are coupled. If the Fermi velocity of the 3d system is 
much larger than in the Id system, the 3<i electrons react practically instantaneously 
to changes in the Id chains. In this case the 3d electrons simply renormalize or "screen" 
the Luttinger parameters in the chains. When Vp ~ v 1 ^ 1 the situation becomes more 
complex, since the dynamics of screening becomes relevant. In this case the charge 
excitations loose their one-dimensional character and become three-dimensional. 

The next step beyond the current work would be to derive an effective model of one- 
dimensional chains coupled via the 3d environment, and to use the methods developed 
and tested in the first chapters to calculate the phase diagram. A major difficulty 
to overcome is the fact that dynamic screening leads to retarded interactions which in 
general cannot be cast into a Hamiltonian formalism, such as the one used in the simple 
model calculations, but rather requires a path integral formulation for the bosonized 
systems. While it is more or less standard to represent a bosonic theory in the path 
integral language, it is not clear how the Klein factors in this formulation can be 
incorporated. 
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Appendix A 

Derivation of the gap equations 



In this appendix we derive the SCHA equations for the Peierls-Hubbard model. The 
derivation for the ionic Hubbard model is similar. For a finite system the Klein factors 
cannot be replaced by their eigenvalues, but we can choose a trial Hamiltonian where 
the bosonic fields are decoupled from the Klein factors as follows 

H tT = Ht + H t f + jjBcsB c b s j (A1) 

with the usual bosonic part, in which the nonlinear terms are replaced by a quadratic 
form 

Ht + Ht = £ I" p \ V -{d x <p a ) 2 + v a9a {d x 9 a f + -^A (A.2) 
^Jo 27r Ua v a g a J 

and a Hamiltonian containing the Klein factors 

h b cs b c b s _ iB cs L(F^F L1 + F^F Ll ) + B c LF^F^F Ll F L ^ 
+B 8 LF+F+F Rl F L ^ + h.c. 

+^{vc9cJc+Vs9sJs)- (A.3) 

Taking the ground state of H tr as a variational state for the full Hamiltonian we obtain 
an upper bond E for the ground state 



E — (Ho)tr + (Hi) tv + (H 2 )tr + (#Peierls)tr, 

E E tr v _R<_r<_« ^ 



dB r dB* dB r ' 



Here £7 tr is the ground state energy of H tT given by Eq. (|A.1|) . and e° r = E® r (B cs , B c , B s )/L 
where E® T (B CS , B c , B s ) is the ground state energy of the Klein Hamiltonian 



B CS B C B S 
tr 
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given by Eq. (|A.3|) . In addition, u = tu/ira, and U = UL/[(2ira) 2 N}. We have used 
the Feynmann-Hellmann theorem (37JI3B] to express the expectation values of the Klein 
factors in terms of the derivatives of e° r , the ground state energy E® T devided by system 
length, with respect to the 5's: 



95 c 

< 
95 s 

< 
dB rx 



(FriFriFliF l1 + h.c), 
(F^F+ F Rl F L ^+ h.c), 
i(F+ F L ^ + F+ l F Ll )+ h.c. 



(A.5) 
(A.6) 
(A.7) 



The parameters A c , A s , B c , B c and B cs are obtained from the condition that E should 
be minimized. For the parameters B c , B s and B cs the minimum condition for E yields 
the following equations: 



uc 



+ 



_ R 

- L/ ( 

Ue-^i) _ B 



J dBl 



u e -m) _ r 



<9 2 e° r 



<9 2 e° r 



0. 



(A.8) 



ue 



B r 



2 o 



9 2 e 



Ue~ 4 ^ - 5, 



J dB c dB, 

|2„0 



+ 



d 2 el 
8B,dB r 



0, 



5, 



2 o 



9 2 e 



dBl 



(A.9) 



we 



5, 



2^0 
tr 



9 2 e 



«9 2 e° r 



+ t/e" 4 ^ 
with the solutions 



5, 



2 o 
tr 



9 2 e 



952 



5 CS 
5 C 
5, 



we 
t/ e -4<^> 



>-(#> 



Minimizing 5 with respect to A c and A s yields 
A? 



2nv c g c 
2nv s g s 



9p° 

_ 4 R ^1 _ p 



-45 



95 L 

< 
'95, 



5, 



9p° 
'dB cs 
del 



tr 



f 95. 



(A.10) 

(A.H) 
(A.12) 

(A.13) 

(A.14) 
(A.15) 
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where we have used that 

d(Al) 2*v a g a [ ■ 1 

These equations are identical to the equations obtained in the thermodynamic limit, 
see Eqs. ()4.28j) - (j4.29j) . if the minimum of e(/c|, k±), given by Eq. (|4.24|) . is replaced with 
e° r , the ground state energy E® r of the Klein Hamiltonian divided by system length. 
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Appendix B 

The trial Hamiltonian in terms of 
bosonic operators 



The bosonic part of the trial Hamiltonian for spinless fermions with static dimerization 
in the thermodynamic limit reads 



o 2tt U 



vg 



After Fourier transformation 



<j>(x) 



it can be rewritten as 



q>0 



V 1 1 a 2 

-q 2 + —A 2 

9 v 9 



-q) + vgq 2 6{q)d{-q) 



If we express 0(±<?) and 9(±q) in terms of the boson operators 



<P{q) 

%) 

'(-?) 





( b Lq ~ b Rq) e 

(bt q + b Rq )e' 



-aq/2 



-aq/2 




(hq ~ b + Rq )e-^ 2 , 
(hq + b + Rq )e-^ 2 , 



2q y ~ 
2q 



(B.l) 



:b.2) 



(B.3) 
(B.4) 
(B.5) 
(B.6) 
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where q > 0, the trial Hamiltonian becomes 



q>0 

E 

q>0 



[bt q b Lg + b + R b Rq ] e 



^g+-A 2 
2tt 2q 



K b t + b t b U e ~ aq + E 

<?>0 



2g 



(B.7) 



This Hamiltonian can be diagonalized with a Bogoliubov transformation of the follow- 
ing form: 



Bi(q) 
B 2 (q) 



ab Rq + /3bl g , 



(3b 



Rq 



ab 



Lqi 



(B.8) 
(B.9) 



with a and (3 parameters to be determined. The operators Bi and B 2 are also bosonic 
operators, i.e. 

[B 1 ,B+] = [B 2 ,B+] = l a 2 -(3 2 = l. 

The operators in the Hamiltonian ()B.7|) are replaced by the new operators B\ and 
B 2 . The restriction that the Hamiltonian in terms of B\ and B 2 has to be diagonal 
yields the second equation for the parameters a and (3 

2a(3{ l4 q 2 + A 2 /2) + ( l2 q 2 + A 2 /2)(a 2 + /3 2 ) = 1, 

with 74 = Up + g i /2irvF and 72 = g 2 /2irvp. The diagonalized Hamiltonian is 



tr 



2 >A> V + A 2 (B+Bx + £+5 2 ) e 

g>0 



aq 



+ E 

<?>o 



- aq 



'7 



-A 2 (l + /? 2 ) + 2/5(/5 7 4 + «7 2 )g 



(B.10) 



The last sum gives a infinite constant which is removed by the normal ordering oper- 
ation. At the end we obtain in the limit a — > 



= £ e(g) (flf-Bi + £+£ 2 ) with e(g) = + A 2 . 

g>0 



(B.11) 



Appendix C 

Quantum theory of Josephson 
junctions 



Here we wish to point out that the study of the Klein factors within the SCHA, for 
the spinless-fermions model, leads to a trial Hamiltonian which is well known from the 
quantum theory of Josephson junctions in the zero-damping limit Considering 
again Eq. (j3.31j) and putting 2pj — > x, we have 

H t B I ^-^l^- 2 +2Bcosx (C.l) 
tr L dx 2 

with periodic boundary conditions for the wave function, ty(x + 2ir) = ^(x). On the 
other hand, the Josephson junction Hamiltonian in standard notation jSHj reads 

H JJ = -4E C -^- - Ej cos <p, (C.2) 



^ + 2tt) = *(<p), (C.3) 

where tp is the order parameter phase difference across the junction, Ej the Josephson 
coupling energy, and E c = e 2 /2C (C is the capacitance) the charging energy. The 
correspondence is obvious. 

Note that the "kinetic energy" in (jC.2|) corresponds to Q 2 /2C, where Q = 2eN is the 
2e times the Cooper pair number, such that N and ip are conjugate variables, [N, <p\ = 
—i. For large junctions, E c Ej (corresponding to large L), the variable ip can be 
treated classically. On the other hand, for small junctions, E c ^> Ej (corresponding to 
small L), the appropriate eigenstates are particle, i.e. Cooper pair, number eigenstates. 
These two limits correspond to two representations of the BCS wave functions, namely 
with a definite phase and with a definite number of Cooper pairs, respectively |3l?] . 
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Appendix D 
Mathieu equation 



The Mathieu equation in its canonical form [3U| reads 

^| + (a-2gcos(2x))y = 0. 
ax z 

This is exactly equation (j3.31j) if we identify 
y(x) — > *(pj), 



(D.l) 



(D.2) 



2LE 



(D.3) 



a 



irvg 



2LB 



(DA) 



irvg 



It can be shown that there exists a set of characteristic values a r (q) which yield even 
periodic solutions and a set of characteristic values b r (q) which yield odd periodic 
solutions, where r = 0, 1, 2, .... If q is real then the sets of characteristic values a r and 
b r have the following properties: 

a) The characteristic values a r and b r are real and distinct, if q ^ 0; a < b\ < a\ < 
b 2 < , q > 0, and a r (q), b r (q) approaches r 2 as q approaches zero. 

b) A solution of the Mathieu equation associated with a r or b r has r zeros in the 
interval < x < tt. 

In Ref. jHlj the asymptotic behavior of the characteristic values a r and b r is derived. 
In particular for q 3> 1 



If we consider the case r = 0, we can estimate the difference E(n) — E(0) needed in 
Eq. dSHSD- We identify E(n) = h(q)q/B and E(0) = a (q)q/B to obtain Eq. (EHnj) . 




(D.5) 




Figure D.l: The lowest eigenvalues of the Mathieu equation as function of q. The a's 
correspond to periodic and the fe's to anti-periodic solutions. With increasing q the 
difference b r — a r _i goes to zero. 
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The Klein Hamiltonian as a 2d tight 
binding model 



In the following we represent the Klein Hamiltonian of the Peierls-Hubbard model, 
h b cs b cBs = iBeaL (F+ Flt + F + p Li ) + BcLF + F + FLiF ^ 
+B S LF+F+ F Rl F^ + h.c. 

+^(v c g c J 2 c +v s g s J*), (E.l) 

in the basis of J c and J s in order to show explicitly that it has the form of a 2c? tight 
binding model for a particle in a harmonic potential. Using the sign convention for the 
Klein factors that we have chosen in section E21 one obtains 



H B CS B C B S = WcsL J2(-l)^\J c ,J s )(J c + 2,J s + 2\ 

+ iB cs L ^(-l)£|j c+ 2, J S )(J C ,J S + 2| 

+ B S L \ J c J s )(Jc Js + 4| + B C L \J C , J S )(J C + 4, J s \ + h.c. 



JciJs J C\J s 



+ ^E J 'I J -^)^^I + ^E J S 2 | J C ,^)(^^I- (E.2) 



4L 4L 

J c-,J s Jc)J s 



For iV c = and iV s = 0, according to Eqs. (|2.56j) and (|2.58|) the quantum numbers J c 
and J s are even and J c ± J s are multiples of 4. In Fig. IE. II the allowed hoppings are 
displayed schematically. 
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The Klein Hamiltonian as a 2d tight binding model 




Figure E.l: Pictural representation of the hopping Hamiltonian (jE.2|) . 



Appendix F 

Analytical solution of the gap 
equations 



In this appendix we present the solutions of the gap equations, both for the Peierls- 
Hubbard model, Eqs. (JPfy - ljOty . and for the ionic Hubbard model, Eqs. (|4.61j) - ([4.62|) . 



Peierls-Hubbard model 

In the case of the Peierls-Hubbard model, for non-zero dimerization the charge and 
spin mean field parameters are given by 



= 4B CS + 8B C = — e-«2>* + % e" 4 ^K (F.l 



2ltV c g c 7T 7T 



4B CS -8B S = — e-<^ - ^ e" 4 ^, (F.2) 



27TV s g s 7T 7T 

where we set a equal to one. If we insert the analytical values for (0 2 )tr and (0 2 ) 
given by Eqs. ()4.32|) - f)4.32|) . we obtain the following set of equations: 

A? Atu ( A„ \ 9c/2 / A, \ 1/2 2U ( A ^ 2£,c 



27TO C # C 7T VAoc/ \Aos/ 7T 



r2 



s/tr, 



(F.3) 
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Analytical solution of the gap equations 



Using Eq. (IF.6I) . we may express y as a function of x 



y = u^X 9 ^ 



2/3 



which we insert in (IF. 51) to obtain 



^4/3 x 2 9c /3 



1/3 



Ux 29c 



We solve this equation in two limits. For u 
constant 

l/(2-2 9c ) 

/ U \ 

X 



U_ 



For u/U C 1 we make the ansatz 
x = x + Cu a , 



(F.7) 
(F.8) 

the solution of Eq. (|F.8|) yields a 

(F.9) 

(F.10) 



with C and a to be determined. Since Cu a is only a small correction we keep it only 
in linear order in the following and obtain 

A c x 2 (l - Uxl 9 ^ 2 ) + 2Cu a (A c x - Ugcxl 9 '- 1 ) 



~4/3 2 5c /3 



2Cg c u c 
x 



1/3 



(F.ll) 

As + UJ 

We consider this equation as a polynomial equation in u which has to be fulfilled for 
every u. Neglecting the term 2Cg c u a /xq <C 1 in the square brackets we obtain the 
following expressions: 



a 



4 
3' 



(F.12) 



C 



(3-4 ffc )/3 

x 



2U(l-g c )(A s + Uy/z 
Finally we reinsert the solution for x into equation (|F.7|) and obtain the result 

A c ( u ) - a c (o) _ f^_y 3 (3_49c)/3 



;f.i3) 



Oc 



X 



A 2 



2U 



-1/3 



g c )U \2nv s g s tt 2 

A S Q) ^ /4 M y /3 /A c (0)V c/3 / A| | 2U 
A 0s \ 7r J \ A 0c / V 27r ^^ 71-2 



-2/3 



;f.i4) 
;f.is) 
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In the other limit, u/U ^> 1 we take the logarithm of Eq. (|F.8J) . 

^4/3 x -4 9c /3 



\n(A c x 



2-2g c 



hi 



U11 + 



{A s + uy/ 3 u'J_ 



and make the ansatz 
x = Bit 13 . 



(F.16) 



(F.17) 



In the limit u ^> U we may neglect the 1 in the logarithm of (jF.lfijl compared to the 
second term and obtain 



B 







1 



3/(6-2 9c ) 



A c (As + uy/ 3 

2 



3 - g c 

Finally the charge and the spin gap are 



:f.i8) 



;f.i9) 



A c (u) 




A 0c 








A 0s 





2/(3- 9c ) 



2/(3-ff c ) 



^Vc9c 

A 2 

^Oc 

A 2 



A 2 , 2 £A _1/(6 ~ 29c) 



2irv s g s 7r 2 / 

2 

n.s 



A 2 , 2 f/\ (4 ~ 9c)/(6_29c) 



2nv s g s n 2 



) 



(F.20) 
.(F.21) 



Ionic Hubbard model 



For the ionic Hubbard model in the case 25 c < B cs the mean field parameters are 
given by the equations 



A' 



2irv c g c 

A 2 
2nv s g s 



ABcs - 8B C = ^ e -^>-^> + ^e- 4 ^), 

7T 7T 2 



;f.22) 

T.23) 



As in the previous case we insert the analytical values of (0^) and (cj) 2 ) and introduce the 
notation A c /A 0c = x, A s /A 0s = y, A 0c /27w c g c = A c , A 0s /2irv s g s = A s , 2A/ti = u 
and 2U/n 2 = U. We obtain 



A c x 2 = Ax 9c/2 y 1/2 + Ux 29c 
A s y 2 = Ax^ 2 y l ' 2 + Uy 2 . 



;f.24) 
;f.25) 
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Analytical solution of the gap equations 



Eqs. (jESl - CT are identical to Eqs. (jF~5j) - (jF~6j) . if u is replaced with A and A s + U 
with A s — U. In the limit A/U ^> 1 we can use the results given by Eqs. (jF.16|) - (jF.21J) 
and obtain 

x2/(3- 9c ) / 2 x3/(6-2 ffc ) , A 2 2C/V l/(6-2 9c ) 

U-^---T , (F-26) 



A C (A) ^ 


(1A 


A 0c ' 




A S (A) ^ 




A 0s 





A 2 Qc / \2ivv s g s TT 2 J 

/2 7rv a \W2(3-9c) / ^ 2 2f/\ (4 ~ 3c)/(6_29c) 



For 2£> c > B cs the gap equations read 

A 2 d2 

— -2- = 8£? c -2-f, (F.28) 

A 2 B 2 
—5- = -8£. + 2-f, (F.29 

which is equivalent to 

2A 2 277 

^ 2 = —rrx^y + — x 29c , (F.30) 

U TT Z 

9A 2 2/7 

A s y 2 = =-z-*y - %y 2 . (F.31) 

U TT 

We assume for x a dependence on A of the form 

x{A) = x + EA' 3 , (F.32) 

with Xo/BA? > 1. If we insert Eqs. (jEHU and (jEH in Eq. (fF~3njl and compare the 
prefactors of same powers of A we obtain for xq, (3 and B the following expressions: 



2U 



A c n 2 



l/(2-2 Sc ) 



x = ( t-^) ) , (F.33) 



/3 = 4, (F.34) 



B ~ A c (l-g c )U 2 (A s + 2U/K 2 ) X ° 29C ~ 1 - (R35) 
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Finally the charge and spin parameter are given by 
A C (A) - A c (0) ^ ,4 ( K , 2U 



2lTV s g s 7T 

2 T /A^y— ) ; 



(l-^ c )f/ 2 A 0c V A 



A,(A) _ A ,/ A^ + 2ay'2/A^0)y»^ 



A 0s V 27ri; a y s 7r 2 / [/ V A 
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